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Abstract. Calculations within QCD (lattice and sum rules) find the lightest glueball 
to be a scalar and with mass in the range of about 1000-1700 MeV. Several 
phenomenological investigations are discussed which aim at the identification of the 
scalar meson nonets of lowest mass and the super-numerous states if any. Results on 
the flavour structure of the light scalars /o(500), /o(980) and /o(1500) are presented; 
the evidence for /o(1370) is scrutinized. A significant surplus of leading clusters of 
neutral charge in gluon jets is found at LEP in comparison with MCs, possibly a 
direct signal for glueball production; further studies with more energetic jets at LHC 
are suggested. As a powerful tool in the identification of the scalar nonets or other 
multiplets, along with signals from glueballs we propose the exploration of symmetry 
relations for decay rates of C = +1 heavy quark states like Xc or Xb- Results from Xc 
decays are discussed; they are not in support of a tetra-quark substructure of /o(980). 
A minimal scenario for scalar quarkonium-glueball spectroscopy is presented. 



1 Introduction 2 

2 QCD expectations for glueballs 4 

2.1 Quantum numbers 4 

2.2 Lattice QCD 5 

2.3 QCD sum rules 8 

2.4 Conclusions from theory 10 

3 Experimental search strategies for glueballs 10 

3.1 Quark-gluon constituent structure of hadrons from their decays 10 

3.2 Enhanced and suppressed glueball production 15 

3.3 Supernumerous states among qq nonets? 16 

4 Spectroscopy of scalar mesons: qq nonets mixing with glueballs 17 

4.1 Route 1: [qq] nonet - glueball mixing around 1500 MeV 17 

4.2 Route 2: light glueball around 1 GeV 18 

4.3 Route 3: two supernumerous states, broad glueball /o(1200 — 1600) ... 20 



2 



5 Properties of /o(500)/(T - comparison with Kq{SOO)/k 21 

5.1 Results on /o (500) /a 21 

5.2 Comparing elastic tttt and Kt: scattering - a problem with the k resonance 25 

6 Intrinsic structure of /o(980) 27 

6.1 Mixing angle from various observables 27 

6.2 Final result on mixing angle 30 

6.3 Possible gluonic contribution 30 

7 Experimental evidence for /o(1370) and /o(1500) - a reassessment 31 

7.1 Production in annihilation 31 

7.2 Search in phase shift analysis of tttt scattering 34 

7.3 Decays of D and B into /o mesons 42 

7.4 Decay of J/i/; 0(lO2O)/o, fo^nn, KK 50 

7.5 Central production in hadron hadron collisions 50 

7.6 Flavour structure of /o(1500) from its decays 54 

7.7 Production in 77 scattering 56 

8 Gluon rich processes in compcirison 57 

8.1 The Pomeron as an example for mixed gluon and quark components . . 58 

8.2 Processes with identified gluons 59 

9 Leading systems in gluon jets 59 

9.1 Rate of neutral clusters at LEP above Monte Carlo expectations 60 

9.2 Proposal for further measurements at the LHC 62 

10 [qq] multiplets and glueballs from symmetry relations: 

first results from Xc decays 65 

10.1 Scalar multiplet according to route 2 65 

10.2 Scalar multiplet according to route 1 - tetraquark expectations 67 

11 Summary and conclusions 67 

11.1 Is there any experimental evidence for the existence of glueballs? 67 

11.2 Arc there supernumerous states in the scalar sector? 68 

11.3 Arc there hints towards a gluonic component in isoscalar mesons? .... 68 

11.4 Future studies of symmetry relations in Xc decays 68 

11.5 A minimal scenario for scalar quarkonium-glueball spectroscopy 69 

List of reprinted figures 79 



1. Introduction 



The existence of "glueballs", bound states of gluons, is a consequence of the 
self-interaction of gluons within Quantumchromodynamics. Such states have been 
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considered to appear in analogy to qq states in a quark-gluon theory of hadrons already 
in 1972 during the early development of QCD [1]. First specific scenarios for glueball 
spectroscopy on the basis of a quark-gluon field theory with hadrons as colour singlets 
have been developed by Fritzsch and Minkowski [2] assuming an analogy between 
massless gluons and photons. Much effort has been devoted during the last 40 years 
into the theoretical analysis as well as the experimental searches of this new type of 
hadrons. 

While a consensus is reached about the existence and some properties of glueballs 
in a world without quarks, the details in the full theory are still controversial. Lattice 
QCD is a formulation of the quark-gluon gauge theory on a space-time lattice and 
it is suitable for describing hadronic phenomena in the non-perturbative regime with 
confined quarks. The theory starts directly from the QCD Lagrangian with no other 
parameters than quark masses and a mass scale. The lightest glueball has quantum 
numbers J^"" = O"*""*" and a mass around 1700 MeV in the theory without quarks but 
could become smaller in the full theory. Alternatively, information on the hadronic 
spectrum has been derived successfully from the "QCD sum rules". In this approach 
some phenomenological input parameters ("condensates") have to be determined from 
experiment but a wide spectrum of quantitative results in hadron physics has been 
established. The predictions from these two theoretical approaches on glueballs will be 
reported and summarized in section 2. 

Most of our attention concerning glueballs is focussed on the scalar sector with 
J^*-^ = O^"*", as the lightest glueball is expected with these quantum numbers. The 
first aim of the experimental search is the identification of scalar nonets of qq bound 
states {q — u,d,s) or other multiplets like those built from "tetraquarks" , according 
to the classification within flavour SU (3) symmetry. Glueballs will then show up as 
super-numerous states. Unfortunately, it is just this scalar sector whose identification 
is still "tentative" while the other multiplets are quite well established [3]. Therefore, 
additional strategics, like the study of "gluon rich" processes have attracted attention. 
We present strategics following such investigations and describe different scenarios for 
a scalar quark-gluon spectroscopy in sections 3 and 4. 

It turns out that the validity of such scenarios depends crucially on the existence of 
certain scalar states and their specific properties. Therefore, in the subsequent sections 
5 - 7 we analyze in greater detail the experimental data on production and decay of 
the isoscalar states which could have a gluonic contribution mixed in and analyze their 
fiavour structure. The Particle Data Group (PDG [4]) has hsted the following isoscalar 
states with J^^ = 0+"^ below a mass of 2 GeV as " estabhshed" : 

/o(500), /o(980), /o(1370), /o(1500), /o(1710). (1) 

Special attention is paid to the reported evidence for /o(1370) which is the cornerstone 
for a particular scenario. This state comes along with /o(1500) nearby in mass, so they 
are discussed together. While attention is paid to data taken already long ago, we are 
witnessing now the advent of high precision results from B factories as well as from the 
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Large Hadron Collider and fixed target experiments at CERN with remarkable future 
potential. 

There are some pros and cons for the different scenarios but finally it is hard to claim 
an ultimate evidence for the existence of a glueball. In the search for a direct evidence 
of gluonic mesons the leading clusters in gluon jets have been scrutinized in several LEP 
experiments. As a result, gluon jets are found different indeed from the expectations of 
Monte Carlo programs without glueballs, otherwise very successful (section 9). 

Finally in section 10, we propose a new tool for identifying flavour multiplets based 
on symmetry relations for decay rates in heavy quarkonium decays, in particular in 
states Xco, Xc2- Violation of such symmetries are indicative of glueball production. First 
available results are promising and should allow for further tests of the existing scenarios. 
The last section contains a summary and our preferred solution for scalar spectroscopy 
with glueball. 

We focus here on questions related to glueball spectroscopy, so there is a large 
variety of experimental and theoretical works on the scalar sector which we do not 
address. There are recent reviews which include complementary material on various 
phenomenological and experimental issues [5] and on other theoretical [6] as well as 
experimental results [7]. 



2. QCD expectations for glueballs 

2.1. Quantum numbers 

In the first specific scenarios for glueball spectroscopy [2] it has been assumed that the 
counting of gluonic states is that obtained for free, massless gluons which have only 
two polarization states, so n gluons can form 2" different states. A simple case is met 
for the two-gluon channel. We consider the colour singlet states which correspond to 
the observable hadronic states. The colour averaged two-gluon systems carry the same 
quantum numbers as the two-photon systems. They have been studied by Landau 
[8] and Yang [9]; the general case is also addressed in a group theoretical analysis by 
Minkowski [10]. The results can be summarized as follows: 

For a single photon there are two helicity states |Aa) with = ±1 corresponding 
to the right and left circular polarization states. For two photons in the cms system we 
consider one photon in +z direction, the other one in —z direction. The states with two 
photons (likewise two gluons in a colour singlet state) have spin J ^ 1 ( "Landau- Yang 
theorem") and C-parity C = +1. Four helicity states |AaAfe) are available from which 
states of definite parity can be formed 

(111) + 1-1-1)); |1-1) ; 1-11); ( |11) - | - 1 - 1) ) 

the first three states have parity P = +1 the last one has P = —1. They form states of 
definite J^^ as can be found by studying the behaviour of these states under rotation 
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and parity P [9] 



0^ 



4++ . . .); 1 2++ I , 4++ . . . ; 3++, 5++ . . . ; |0-+| , (2"+, 4"+ . . . (2) 



Here the states with J^^ — 2++ and descendants correspond to the states Jz — ±2 
above. Those states which are formed by "constituents" in a relative S wave are expected 
to be lowest in mass and are enclosed by a box in (2); the higher spin excitations 
correspond to even orbital wave functions. 

Three gluon systems can have both C = +1 and C = —1. For the ground states 
with relative S waves one finds [2] 

C = +1 : J^^ = 1++ 

C= -1 : J^^ = 1—, 3—. 

2.2. Lattice QCD. 

Because of "asymptotic freedom" QCD interactions become weak for short distances 
where they can be analyzed within perturbation theory at small coupling between quarks 
and gluons. This method fails at large distances where the interactions become strong. 
Lattice QCD is a formulation of the theory on a space-time lattice and it is suitable 
for describing hadronic phenomena in the non-perturbative regime [11] (see [12] for an 
introduction and [13] for a recent review). Physical quantities can be computed from 
the QCD gauge and fermion actions Sq and Sp with the quark masses niq as parameters 
and a suitable mass scale. 

The theory can be defined on a Euchdean space-time lattice with fermions on 
the lattice sites and gauge fields connecting the sites. The lattice spacing a acts as 
an ultraviolet cutoff which provides a gauge invariant regularization of QCD. Physical 
observables can be obtained from an average over the relevant configurations of quark 
and gluon fields according to the Boltzmann weight e~^''~^° within a space time volume 
L^T and the physical results are obtained ultimately by extrapolation to the continuum 
limit a — )> 0. 

Hadron masses are calculated from 2-point correlation functions Cij{t) for operators 
Oi relevant for the hadrons considered 

Cij{t) = ^JdijjdijJ dUe-^''-^^{0\Oi(tyOj{0)\0) (3) 

in terms of a path integral over fermionic and gauge field variables ijj, ijj and U with 
normalization Z. The correlation functions decrease exponentially for large Euclidean 
times t 

Cij{t) = cfcf e-^-* (4) 

where the term with the lowest energy Eq dominates and determines the mass of the 
ground state. All masses are obtained in units of lattice spacing a so that only ratios of 
masses are predicted. In order to relate to absolute scales a suitable physical mass scale 
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Figure 1. Mass spectrum of glueballs (in GeV on r.li.s.) for different quantum 
numbers PC according to the quenclied lattice calculations (figure from [18]). 



has to be taken as input, such as the "string tension" or the "Sommer scale" l/rg ~ 400 
MeV [14]. 

As examples of recent lattice calculations from first principles, we mention the 
results in full QCD on the conventional light hadron spectrum by the Budapest- 
Marseilles- Wuppertal Collaboration [15] who has calculated the masses of the baryon 
octet and decuplet states as well as the masses of some light mesons within a few 
percent of accuracy. Here the masses of tt, K and S particles have been used to fix the 
masses of light and strange quarks at their physical values as well as the overall mass 
scale. Another result, obtained by the "Hadron Spectrum Collaboration" [16] concerns 
the spectrum of lightest and the first excited isoscalar meson states which includes 
quark-annihilation contributions. Remarkably, the mixing pattern of these mesons is 
reproduced close to observations. 

More difficult to compute is the spectrum of glueballs in full QCD, as these states are 
heavier and therefore need higher statistics, in particular the scalar states with vacuum 
quantum numbers have extra contributions difficult to disentangle. In full QCD there 
is a mixing of gluonic and fermionic degrees of freedom, correspondingly one inserts 
gluonic and fermionic operators for the relevant correlation functions. For sufficiently 
light quark (pion) masses the glueball can decay into a meson pair which has to be 
included in the consideration as well. 

The spectrum of glueballs has been calculated at first within the pure (Yang-Mills) 
gluon theory without quarks ("quenched approximation"). The lightest glueballs are 
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found for J^^ = O"*""*", 2++ and O""*" which correspond to the S wave ground states of the 
three different spin configurations hsted in (2). Masses have been obtained by several 
groups in good agreement within around 10% by Bah et al. [17], by Morningstar and 
Peardon [18] and with a similar calculation but for larger lattices and volumes by Chen 
et al. [19] at a lattice spacing of around a = 0.1 — 0.2 fm (see Table 1). The results by 
Morningstar and Peardon [18] for states below ~ 5 GeV are shown in figure 1. The three 
S wave states above are indeed the states of lowest mass, the hghtest states formed by 
three gluons with J^^ — l+~ and 1 for C = — 1 are found at around 3 GeV or higher. 

Table 1. Glueball masses (in McV) in quenched lattice approximation; an additional 
error of ~ 5 % should be added for the scale error of 1/ro = 410(20) MeV [18]. 



jPC 


Bah et al. [17] 


Morningstar et al. [18] 


Chen et al. [19] 


0++ 


1550 (50) 


1730 (50) 


1709 (49) 


2++ 


2270 (100) 


2400 (25) 


2388 (23) 


0-+ 


2330 (260) 


2590 (40) 


2557 (25) 



The infiuence of dynamical qq contributions in full QCD on the mass of flavour 
singlet gluonic mesons has been studied by Hart and Teper (UKQCD collaboration 
[20]). In this exploratory investigation the lattice spacing was around a ~ 0.1 fm and 
the light quark masses about one half the strange quark mass; an extrapolation to 
physical limits was not possible. A significant suppression of the scalar glueball mass by 
the amount ~ 0.84±0.03 with respect to the quenched value was found. With increasing 
quark mass this suppression should disappear and pure gluo-dynamics be restored. As 
this effect was not observed, not even for quark masses twice the strange quark mass, no 
real significance has been attributed to this suppression effect. Furthermore, the spin 2 
glueball did not show any substantial deviation from the quenched result. 

A different conclusion has been reached by Hart et al. (UKQCD [21]) within a 
study using both gluonic and fermionic operators to create the flavour singlet states. 
There have been Nf — 2 degenerate sea-quarks using a sample of a total of ~1500 
conflgurations. At comparable lattice spacing the mass of the lightest gluonic meson is 
reduced considerably with respect to the unquenched result from 1600 McV to about 
1000 MeV. The results on the two lightest 0"'""'" mesons are interpreted in terms of a 
maximal mixing of the qq and gg states into the physical flavour singlet mesons around 
1000 and 1600 MeV. 

More recently, a study with higher statistics and including strange quarks has been 
presented by Richards et al. (UKQCD [22]). They use measurements at pion masses 
of 280 and 360 MeV and a spacing of a = 0.123 fm and a = 0.092 fm with 3000 and 
5000 configurations respectively. In their analysis gluonic and tttt operators have been 
included but no fermionic operators. This approach did not show any essential difference 
between unquenched and quenched results for the lightest O'^"'" state. A mixing of the 
glueball with the qq state as found in [21] has not been considered explicitly. 
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Similarly, for the 2++ and O""*" states the results on the masses show little difference 
between quenched and unquenchcd calculations and the same conclusion has been drawn 
by the same collaboration (UKQCD) for glueballs of yet higher mass [23]. A summary 
of the unquenched calculations is presented in Table 2. It should be noted that the 
scale parameter vq is different in quenched and unquenched calculations, but the values 
actually used in [18] and [22] differ by < 3%. 



Table 2. Glueball masses from lattice QCD: listed are the ratios of masses from 
unquenched over quenched calculations. 



jPC 


Hart and Teper [20] 


Hart et al. [21] 


Richards et al. [22] 


0++ 


0.84 (3) 


0.63 


1.03 (3) 


2++ 


1.03 (3) 




0.98 (10) 


0-+ 






0.97 (2) 



For the future studies a yet higher accuracy of the simulations is suggested in [22] , 
also smaller quark masses should be used to allow for a more realistic decay of the 
glueball into two pions. Furthermore it is argued that a definitive calculation required 
a continuum extrapolation, and the inclusion of fermionic operators [23]. In that way 
the mixing problem between gg, qq and tttt states could ultimately be clarified. 



2.3. QCD sum rules 

An alternative approach to obtain information on the hadronic spectrum within QCD 
and on the glueball properties in particular is based on the "QCD sum rules" [24] 
(see [25] for a pedagogical presentation and [26] for a recent review). It starts from 
a correlation function in Minkowski space involving gluonic or fermionic operators Oi 
according to the considered J^'-^ quantum numbers which is evaluated in the deeply 
space hke region 

U{Q') =tj d'x e"^^ (0|T{Oi(a;)O2(0)}|0) {Q' ^ -q' » A^^^). (5) 

In an application of the operator product expansion the short distance contribution to 
the correlation function is calculated from QCD perturbation theory whereas the long 
distance contribution is given in terms of "vacuum condensates". These parameters, 
ultimately, have to be determined from experiment, of special importance are the 
condensates of lowest dimension, the quark and gluon condensates (qq) and {^GIvG'^^^) . 
One can estimate the low energy hadronic spectrum by matching the correlator (5), as 
given by this expansion, with a dispersion relation of the type 

in general, (6) has to be modified by subtractions depending on the high energy 
behaviour of the absorptive part p{s) = Imn(s) which adds terms with powers of 
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in (6). The absorptive part can be represented as a sum over the low lying hadronic 
states and a continuum contribution, in the simplest case of one resonance as 

p{s) = p^^'is) + 9(3 - So) Imn«^^(s), -/^^ = fMs - mlJ, (7) 

TT 

where sq denotes the threshold of the QCD continuum. For better sensitivity to low Q"^ 
of (6) one performs a "Borel transformation" which calculates the integral 

L^''){r) = - / dtt''e-nmSl{t) (8) 

with index A; = — 1,0,1,.... In this way a set of sum rules for index k is obtained. 

The successful application of this formalism to a large variety of hadronic 
phenomena concerning conventional mesons and baryons has given a good confidence 
to this approach besides providing the necessary phenomenological parameters. 
Consequently, this scheme has been applied to the analysis of glueballs as well. 

It has been suggested that the valence gluons bound in a glueball couple much 
stronger to the vacuum fields, than, for example, the quarks in the p meson where a 
description in terms of mean vacuum fields is sufficient [27]. Ultimately, this strong 
coupling to vacuum fields leads to the higher mass scales for glueballs. Particularly 
strong effects from vacuum fields occur in spin zero systems of low energy, a problem 
which has been approached by the inclusion of instanton field configurations. Along this 
fine of thought the radius of the 0++ glueball is only 0.2 fm to be compared with 0.5 fm 
for the p meson [28, 25]. 

The analysis of sum rules with lowest indices k = —1, (referred to as "subtracted" 
[27] and "unsubtracted" [29] sum rules) by Narison and Veneziano [30] leads to a 
consistent solution including two gluonic scalar states with masses (see [31] for an 
update) 

Mgbi ~ 0.9 - 1.1 GeV, Mgb2 ~ 1.5 - 1.6 GeV (9) 

Furthermore, for the state around 1 GeV, the application of certain low energy sum 
rules predict a large width into tttt with [30, 31] 

r(M,M^,+,-) ~ 0.7 GeV. (10) 

In an extension of this result a phenomenological scheme has been suggested for 
/o(500 — 1000) /o" I, /o(980) and a sequence of isoscalar mesons heavier than 1 GeV 
realizing a mixing of qq and gluonic states as in (9) and a radial gluonic excitation [32] . 

An approach emphasizing the role of instantons by Forkel [33] suggests a solution of 
sum rules with a single scalar state at a mass of Mg6(0++) ~ 1.25 ± 0.2 GeV, in between 
the previous results (9). The same approach has also been applied to the pseudoscalar 
gluonium and gave the result Mgi,{0~~^) — 2.2±0.2 GeV, similar to the lattice predictions 
mentioned above. 

I The higher mass 1000 MeV refers to the "Breit-Wigner" or "on-shell" mass which are more 
appropriate for sum rules (see section 5.1.1). 
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More recently, a sum rule study of the mixed scalar system of gluonium (gg) and 
quarkonium (qq) has been presented by Harnett et al. [34] which also included instanton 
contributions. The analysis of the diagonal qq and gg as well as the non-diagonal qg sum 
rules yields a consistent solution when using two mass states in (7) at approximately 
1 GeV and 1.4 GeV, similar to (9). In this solution the couplings of the two hadronic 
states to qq and gg "constituents" is comparable, corresponding to a near maximal 
mixing, with a slight preference of the heavier state to couple to gg. This would imply 
that there is only one glueball and one quark- antiquark state in this scheme for the two 
hadronic states. The heavier mass state could be related to /o(1500) and the lower mass 
state to either /o(500 — 1000)/(T or /o(980), both higher and lower mass states being 
mixed of gg and qq. 

2.4- Conclusions from theory 

There is a general agreement in that the lightest gluonic state has quantum numbers 
J^*^ = O"*""*". One state is located around 1.4-1.7 GeV. A second state is suggested in 
sum rules and in certain unquenched lattice calculations near 1 GeV. In these analyses 
there is a strong, almost maximal mixing of the hadronic states near 1 and 1.6 GeV 
into the qq and gg components. This agreement between both approaches is striking 
but not universally accepted. In particular, the most recent lattice calculation does not 
find this low mass state, but, given the computational limitations, its existence cannot 
be excluded either. A definitive calculation requires a higher statistics, a continuum 
extrapolation and the inclusion of fermionic operators. 

The next heavier states are expected with quantum numbers J^'^' = ^,2"*"+ and 
with masses > 2 GeV. States with other quantum numbers including exotic, i.e. non-gg 
ones, are expected for higher masses. Experimental analyses have been difficult so far 
in this mass region. Therefore, the search for the scalar gluonic states looks particularly 
promising despite the experimental and theoretical uncertainties. 

3. Experimental search strategies for glueballs 

In the following we will mainly focus on the scalar glueball which is expected to be 
the lightest one and which has been addressed in most research studies. The isoscalar 
states with J^*^ = listed as "established" by the PDG [4] below a mass of 2 GeV 
are already listed in (1). All these states have been considered as glueball candidates or 
mixed glueball/quarkonium states in some of the models or theories. First we discuss 
some conventional strategies for distinguishing glueballs. 

3.1. Quark-gluon constituent structure of hadrons from their decays 

The 2-body strong decays of an unstable hadron reflect its constituent structure. Under 
SU (3) fi symmetry of the strong interactions a light quarkonic meson M{qq') decays by 
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creation of quark anti-quark pairs from the vacuum with equal amphtudes 

|0) ^ \uu) + \dd) + \ss). (11) 

into a superposition of meson states M{qqi)M{qiq'). A glueball can decay by creating 
qq pairs and form quarkonic mesons or by gluons and form gluonic mesons. An expUcit 
model for such decays has been proposed in [35] and a similar approach in [36, 37] with 
small differences for glueballs. The role of OZI suppressed processes has been studied 
by Zhao [38]. 

The flavour symmetry can be broken in different ways. 

• Strange quark suppression: because of the heavier strange quark mass the ratio 

p= {ss\V\0)/{dd\V\0) (12) 

of matrix elements for the creation of ss versus uu or dd is suppressed, i.e. p < 1. 
Empirically, p > 0.8 for estabhshed nonets, for tensor mesons p — 0.96 ± 0.04 is 
found in [35] using form factors, in [36] p — 0.7 — 0.9 without form factors. 

• Form factors: they may provide momentum dependent effects, representing, 
for example, inelastic channels at higher momenta or the influence of angular 
momentum barriers. In [35] the choice is 

|F,,(g-)p = exp(-gV8/32) (13) 

which multiplies the phase space and threshold factor 

SAq) = q''^' (14) 

for angular momentum i and momentum ? = |^ with ^ — 0.4 GeV. 

• Chiral effects: For the decay of a scalar glueball G the amplitude M{G qq) oc ruq 

(quark mass m^) according to a perturbative analysis using chiral symmetry which 
involves the intermediate process gg — )■ qq [39]. In consequence, the decay branching 
ratio rK-K = B(G K~) / B{G — ?■ tt+tt") is enhanced over rxn = 1 for flavour 

symmetry, whereas strange quark suppression would lead to rx-K < 1- 

3.1.1. Quarkonium (qq) decays. In our comparison with glueball decays we consider 
here mainly isoscalar quarkonium states \qq) which are generally mixed from non-strange 
\nn) — -^\uu-\- dd) and strange quark \ss) components 

\qq) = cosQ;|nn) — sinQ;|ss). (15) 

Alternatively, the mixing can be described in terms of the SU{3) fi singlet and octet 
components |1) — -^luu + dd + ss) and |8) = -^\uu + dd — 2uu) . The mixing angle a in 
the quark-flavour basis is related to the nonet mixing angle 9 in the singlet-octet basis 

by 

a = e + arctan ^2 ~ ^ + 54.7°. (16) 
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For a = (7r/2) we obtain pure nn (ss) states ("ideal mixing") as approximately found 
for the vector mesons p, ui^cf). The mixing of pseudoscalar mesons is written as 

I77) — cos (f)ps\nn) — sva.(f)ps\ss) , (17a) 
I77') = sva.(t)ps\nn) +cos0ps|ss). (17&) 

According to a recent review [40] , the different analyses from radiative decays of vector 
and pseudoscalar mesons yield values 42° {Ops — —13°), not far away from the pure 

flavour singlet-octet states {9 — 0, — 7r/2 or = 54.7°, — 35.3°J). On the other hand, 
the possibility of a gluonic component of 77' remains ambivalent. Such a component can 
be added to 77', for example, by introducing an additional mixing angle 0Gps 

1 77') = cos (pGps {sm.(pps\nn) + cos 4)ps\ss)) + sin (pcps \gg)- (18) 

For the last term with Z^/ = sin^^p^ one finds typically \Z^i\'^ 0.1 but larger values 
up to l^^'P < 0.5 are reported as well [40]. A more general approach with two mixing 
angles describing the mixing of meson states at the hadronic scale and also the mixing 
of decay constants at short distances is discussed in [41, 42]. Our analysis of Xc decays 
in the next subsection suggests that the conventional state mixing approach as above is 
appropriate in that case. 

According to the model [35] the partial width Fjj of a quarkonium state into a pair 
of mesons MiMj is computed from 

^ X \F^jm' X SM (19) 

with (13), (14); it depends on the above parameters a and p through the couplings 

7j = Q,-|M,,f (20) 

where Mij is the decay amplitude and Cij is a weighting factor for the different charge 
states. For an isoscalar (qq) decay these couphngs 7|^- are given in table 3 for — 
[35, 36]; the dependence of jfj on a is shown in figure 2. 

3.1.2. Glueball decays. Bound gg states ("gluonia") decay through a two gluon 
intermediate state where several mechanisms can be envisaged [35] (see figure 3) 



a) gb^gg^qqq'q' M M' (21a) 

b) gh^ gg^ gggg r]{r]') r]{r]') or gb gb (21&) 

c) gb ^ gg ^ qq ^ qqq'q' M M' (21c) 

In the decay rates of these processes the ratio of amphtudes 

R = {s-s\V\g)/{dd\V\g) (22) 
I Previously, a favoured value was quoted as (j)ps = 37.4°, 6ps = —17.3° [3]. 
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Figure 2. Couplings for isoscalar (qq) meson decays as function of mixing angle 
a as in (15) for SU{3)fi symmetry (p — 1) from table 3 with tq — for (j)ps — 37.4° 
(figure from [35]). 



R 




•;,;--;.<c:5 n(n') 

G <"> 
•:■-■■><==> nW) 

R 




Figure 3. Processes contributing to gluonium decay (21a-21c) (figure from [35]). 



Table 3. Couplings 7^^ for decay of an / = meson /o with mixed (qq) and (gg) 
components [35, 36]; fiavour mixing angles a as in (15), (p = (j)ps as in (17a); glueball 
decay according to (21a) [35]; 5/7(3)/; breaking parameters p,R as defined in (12) 
and (22); weights Cij included in jf^; decay couphngs gfj — g'^jfj/4: for (qq) and 
ra = ^3/(2 + i?2)G/g for (gg) with g, G from (24) (based on [36, 37]). 



decay channel couplings 7?- for |/o) = cos^gI??) + sin(f>a\gg) weights 
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- i?2 sin2 0))2 


1 
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-i?2))2 
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(cos a SVC? (j) — \/2p sin a cos^ + re (sin^ (j) + 


i?2 C0S2 4))f 


1 
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appears for the qq production from a gluon, where R"^ — 1 for SU (3) fi symmetry. For 
the decay of a mixed state 

l/o) = cos (f)G\qq) + sin (f)G\gg) (23) 

with \qq) mixed as in (15) and \gg) decaying as in (21a) into pairs of pseudoscalars we 
write the decay couphngs as [37] 

9ij-9 d!^ + G d^f; g ^gocos cj>G, G = Go sin 0G. (24) 

The couphngs ^fj for the decay of the mixed state with the factor g'^/A taken out and 
with ro = V3/(2 + R^) G/g are given in table 3. For the glueball decay with flavour 
symmetry [R — 1) one finds from the term with ra oc G the reduced partial widths ^yfj 

-f^jigb ^TTTT : KK : r]r] : riri' -.7]' r]') = 3:4:1:0:1 (25) 

according to the statistical weights. 

The two-gluon intermediate process also dominates the decays of scalar heavy 
quark, cc, bound states with C parity C = +1 like Xc in the perturbative analysis and 
one expects decay rates as in (25). An update of this test for spin J — 0,2 charmonia 
XcO,Xc2 [35] is carried out in table 4. We list the measured branching ratios B into 
pairs of pseudoscalars according to the PDG [4] , where recent contributions come from 
CLEO [43] and BES [44]. From the model (but without dynamical formfactor Fij{q^), 
i.e. /3 — )■ oo) one expects B oc In table 4 we also show B after division 

by cg^^^^, and after normalization to that value for the tttt decay. These data are to 
be compared with the model predictions for |A'/p = in table 3 for ((?(?), i.e. for 
g/G ^ 0, where recall g and G are given by (24). These expectations are given in table 
4 as well for different values of the strange quark parameter R in (22). We observe, that 
the expectations from SU{3)fi symmetry (i? = 1) are fulfilled for all decay processes 
of Xc0,2 within 10 - 20% (only the upper bound for Xc2 — > v'v' is a bit low). So there 
is no strong need for the additional gluonic processes (21&) affecting r],ri'; nor is there 
any strong need for dynamical form factors Fij[q^)] the weakness of angular momentum 
barrier factors (of Blatt-Weisskopf type) may be due to the small radius of charmonium 
states of Oil/rric). 

Looking at the substructures at the 10% level we note that for the tensor state 
Xc2 the observed values for B/{cq^) would prefer a slightly reduced R > 0.9. In case 
of the decays of the scalar Xco values of B/{cq) for kaons are slightly enhanced by 10- 
20%. This effect could be a remote consequence of the chiral enhancement in scalar 
decays discussed above from the contribution of the intermediate process gg ^ qq as 
considered for scalar glueballs in (21c) (the argument is also applicable for a decaying 
scalar cc state). The effect can be simulated by choosing R> 1 and we give in the last 
line of table 4 the predictions for R — 1.1 which would fit the data for Xco generally 
better. 

A modification of the SU (3) fi symmetry for decay rates (with i? = 1 in table 4) is 
expected in schemes which include state mixing at hadronic scales and mixing of decay 
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Table 4. Branching ratios B of Xco(3415) (0++) and Xc2(3556) (2++) into pairs of 
pseudoscalars; also B/{cq) or B/{cq^) resp. with charge weight c and momentum q 
after normahzation to the value for tttt (data as listed in [4]), together with the resp. 
theoretical expectations |Mp = 7|,/c for process gg in table 3 for different strange 
quark parameters R (mixing angle (f)ps ~ 37.4° for r], r}'). 
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3.14±0.18 


3.03 ±0.21 


< 0.24 


2.02 ±0.22 
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R=0.95 


1. 


0.9 


0.9 


0.93 


0.002 


0.88 


R=0.90 


1. 


0.81 


0.81 


0.86 


0.008 


0.77 


R=l.l 


1. 


1.1 


1.1 


1.08 


0.002 


1.13 



constants at short distances [42]; in that case the ratios of the Xc decay rates deviate 
considerably from unity, for example B{xco ~^ VV)/B{Xco ~^ 7r°7r°) = 1.9 and the same 
follows for the ratio {r]'r]') / {n^n^). The data appear to be closer to the conventional 
meson state mixing approach as seen in table 4. 

The remarkable success of the simple scheme based on SU{3)fi symmetry with 
phase space correction for decays into the same nonets suggests further tests with other 
nonets in the search for extra gluonic states, as will be discussed in section 10. 

3.2. Enhanced and suppressed gluehall production 

Not only the decay, also the production properties are characteristic for the intrinsic 
structure of a hadronic state [45, 46]. 

3.2.1. Gluon-rich processes. There are processes which provide a "gluon rich" 
environment with enhanced probability for glueball production, see figure 4 for examples. 
These processes have been extensively explored experimentally in the past [5, 7]. We 
give here a survey first over different processes and come back to specific problems later: 

(i) Radiative J/tp or T decay: In the perturbative approach the (cc) or other heavy 
quarkonium J^^ — 1 bound states decay predominantly through {QQ) 3g. 
Alternatively, the radiative decay (QQ) — >■ 7 + {2g) is possible with formation of 
an intermediate gluonium (gg) state {{QQ) — > 7 + gb). 



16 




Figure 4. Processes favouring glueball G production, where J/f/' and T are 
respectively the lowest mass cc and bh mesons with J^'-^ = 1~~ . 



(ii) Central production of mesons: In double diffractive high energy processes the 
incoming hadrons scatter with small momentum transfers and carry on the initial 
valence quarks. In Regge theory this process is dominated by "double Pomeron 
exchange". If the Pomeron is viewed as a gluon dominated object then glueball 
production is enhanced in this reaction {pp p gb p). The different contributing 
processes within QCD have been discussed in [47]. 

(iii) pp annihilation: The annihilation of quarks may proceed through intermediate 
gluons and the formation of glueballs (pp gh + M). 

(iv) Decay of excited heavy quarkonium F*^") to ground state Y : In the example 

-^Y -\- X the hadrons X are emitted from intermediate gluons and therefore 
could be formed through an intermediate glueball. 

(v) Decay of heavy quark b — > sg: This QCD process (through "penguin" diagram) 
may hadronize involving a glueball according to S — > Kgb [48] . 

(vi) Leading particle in gluon jet: In analogy to the fragmentation of the primary quark q 
of a g-jet into an energetic meson M{qq') which carries q as valence quark, there may 
be the fragmentation of the primary gluon of a gluon jet into an energetic meson 
M{gg) which carries the initial gluon as valence gluon g ^ gb + X (section 9). 

3.2.2. Suppression of glueballs in 77 processes. Having neutral constituents a glueball 
couples to photons only through loop processes and then it is suppressed in 77 reactions. 

3.3. Supernumerous states among qq nonets? 

Mesons with light quark constituents (-u, rf, s quarks) are classified in nonets of 3 x 3 
states (octet+singlet). A well known example is the pscudoscalar nonet of lowest mass 
with TT, i^T, T] near flavour octet and rj' near singlet. It is the aim of meson spectroscopy 
to estabhsh the appropriate classification of mesons. They should fit into nonets of qq 
states - possibly, there are also exotic states like tetra-quark qqqq or hybrid qqg states. 
If there are glueballs in addition there should be supernumerous states which do not fit 
into a nonet classification of the meson spectrum. 



17 



4. Spectroscopy of scalar mesons: qq nonets mixing with glueballs 

In the scalar sector different scliemes liave been proposed for a spectroscopy witli 
glueballs for the mesons below 2 GeV including the isoscalars of (1). 

Table 5. Spectroscopy with glueballs, route 1: glueball near 1500 MeV. 



/ = 




1=1 




ii:*(1950) 


/o(1710) 
/o(1500) 
/o(1370)? 


ii'*(1430) 


ao(1450) 


3 isoscalars: qq nonet 
+ glueball M~ 1.6 GeV 


/o(980) 

/o(500)/a 


if*(900)/K? 


ao(980) 


light nonet 
qq or qq - qq 
or KK 



4.1. Route 1: {qq) nonet - glueball mixing around 1500 MeV 

Lattice QCD in quenched approximation predicts the lightest glueball in the mass range 
around 1600 MeV (see section 2.2). The discovery of the rather narrow /o(1500) in pp 
annihilation by the Crystal Barrel Collaboration [49] and by GAMS [50] suggested 
at first the possibility of a gluonic state. A closer inspection of decay mechanisms 
and branching ratios led Amsler and Close [35] to propose a mixing scheme including 
/o(1370) and /o(1710) with a gluonic component as well §. In table 5 we list the meson 
states ordered according to mass and isospin. Between 1.0 and 1.9 GeV we find the 
states fitting into a qq nonet but with one supernumerous isoscalar which suggests the 
presence of a glueball. 

With more experimental results available and with further development of models 
various mixing schemes for the three isoscalars have been proposed (for a further 
discussion, see [7]). As an example, we present the result [51] which includes data 
for 77 processes. The physical mass states foirrii) are decomposed into the qq and gg 
states or, alternatively, into the 5'[/(3) eigenstates |1) = {\uu) + \dd) + \ss))/\/?> and 
|8) = {\uu) + \dd) — 2|ss))/-\/6 and gluonium \gg) 



l/o(1370)) 
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-0.50\ 


\nn) 
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-0.50\ 


|1) 


|/o(1500)) = 
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-0.61 


0.61 


■ \^^) — 


0.00 


0.75 


0.61 


• |8) 


l/o(1710)) 


^0.22 


0.76 


0.60 y 


\99) 


^0.62 


-0.49 


0.60 y 


\99) 



One observes that in this scheme /o(1500) is a mixture only of octet qq and glueball; the 
glueball component is distributed over all three /o states (last column in mixing matrix) 
with about equal amounts. In alternative schemes either /o(1500) can be dominantly 
gluonic [52] or /o(1710) while /o(1370) is near singlet [53, 52]. 

§ In the scheme [35] the existence of the heavy state (now /o(1710)) actually has been predicted. 
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A possible problem with this mixing scheme are the doubts related to the very 
existence of /o(1370). This resonance is included in global fits of parametric model 
amplitudes to a variety of channels but it is not seen in any model independent bin-by- 
bin phase shift analysis. The PDG under this entry does not provide any established 
result on branching ratios or ratios thereof - contrary to the well established nearby 
/o(1500). We come back to this problem in more detail in section 7. 

The light scalar mesons below 1 GeV can be grouped into another nonet (see 
table 5) which includes the broad states /o(500)/(T and K*{SOO)/k,, besides the narrow 
/o(980) and ao(980). This meson nonet can been built from qq as usual, see, for example 
[54, 55, 56]. Alternatively, it can be given a substructure of diquarks {qq — qq) according 
to the scheme by Jaffe [57] within the MIT bag model where the nonet is constructed 
as the direct product of the anti-triplet of diquarks {ud, us, ds) and the triplet of anti- 
diquarks {ud, us, ds). This construction explains naturally the ordering in mass of the 
states where Gq, Jq have two, k, has one and a has no s{s) quark constituent 

a+ = [su] [sd], /o/a° = {[su] [su] ± [sd] [sd])/V2, a~ = [sd] [su] 

k'^ — [su] [ud], = [sd] [ud], = [ud] [sd], kT = [ud] [su] (27) 

(7 = [ud] [ud] 

Various applications of the 4-quark model for the hght scalars have been considered for 
decay and production processes, in particular with photons [58, 59]. A picture of these 
scalar mesons as a mixture of tetra-quark states (dominating in the light mesons) and 
heavy (qq) states (dominating the heavier mesons) has been proposed in [60] based on an 
instanton induced effective Lagrangian theory. The mixing between the mass eigenstates 
/o and cr and the ideally mixed states in (27) is typically assumed small to be consistent 
with the near mass degeneracy of /o and Oq. For example, the corresponding mixing 
angle is < 5° in [60]. 

The hght scalar mesons can be estabhshed by the motion of phase shifts through 90° 
in the appropriate two-body processes, except for the k with a phase in i^Tr-scattering 
staying below 40° (see section 5.2). 

4-2. Route 2: light glueball around 1 GeV 

The QCD results discussed above do not allow a definitive conclusion about the mass 
and properties of the lightest scalar glueball. A light glueball near 1 GeV is obtained 
from QCD sum rules; based on low energy theorems also a large width is suggested 
which points to /o(500), sometimes also referred to as /o(500 — 1000) according to the 
heavier Breit-Wigner mass of ~ 1000 MeV (see below and section 5.1.1). A light glueball 
near 1000 MeV is also expected from certain lattice calculations [21] (see section 5.1.1), 
both states /o(500 — 1000) or /o(980) are possible candidates. 

In view of the theoretical ambiguities at the end of the 90's, a largely 
phenomenological approach to light scalars including a glueball has been pursued by 
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Minkowski and Ochs [61] without reference to the above specific QCD results. General 
properties, especially mixing patterns are studied basing on an effective action for sigma 
variables [62]. The lightest qq nonet is constructed from the list of well established 
resonant states, see table 6. In particular, the existence of /o(1370) as an extra state has 
not been accepted: the peaks appearing in the tttt mass spectrum have been interpreted 
as caused by a broad object centered around 1 GeV interfering destructively with the 
narrow resonances /o(980) and /o(1500) ("red dragon" phenomenon). This broad object 
has been related to /o(500): after separation of /o(980) the tttt elastic phase shifts pass 
90° near 1000 MeV and can be represented locally by a "Breit-Wigner" resonance. A 
typical fit for '7o(1000)" yields (see also [63]) 

Mbw ~ 1000 MeV, Vbw ~ 700 MeV. (28) 

The determination of the amplitude pole from the phase shifts is non-trivial and we 
come back to this question in section 5.1.1. 

The hghtest qq nonet in that scheme includes the isoscalars /o(980) and /o(1500) 
with a mixing similar to rj' — r). This has been motivated by the observed decay rates 

J/ip ^ uj,(j) + X which project out the strange and non-strange quark components of 
X|| and are found comparable for X — rj' and X — /o(980). According to PDG [4] (in 
units of 10"^): 

B{J/^l) ^UT]') = 1.82 ±0.21; B{J/^ ^ (pr]') = A.O ±0.7 (29) 
S(J/V' ^ w/o(980)) = 1.4±0.5; S( J/^ ^ (/)/o(980)) = 3.2 ± 0.9 (30) 
We write the flavour mixing of these scalar mesons as 

|/o(980)) =sm(j)sc\nn) + cos (j)sc\ss), (31) 
|/o(1500)) = cos (f)sc\nn) - sin^^dss), (32) 

or inversely 

\nn) = sin(/),e|/o(980)) + cos (/),e|/o( 1500)) (33) 
\ss) = cos(/.,e|/o(980)) - sin(/.,e|/o(1500)). (34) 

II It is assumed that the "singly disconnected diagram", where the qq is produced in a single loop 
through one intermediate gluonic exchange, dominates over the "doubly disconnected diagrams" [40] . 



Table 6. Spectroscopy with glueballs, route 2: light glueball around 1 GeV. 
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For estimates, the flavour composition \uu, dd, ss) has been approximated by 



7j', /o(980) ^ 1 1, 1, 2)/VQ; v, /o(1500) ^ | 1, 1, -1)/V3. 



(35) 



The octet formed from ao(980), K*{U30) and = /o(1500) is found to follow the 
Gell-Mann-Okubo mass relation (using masses 984.7, 1412 and 1507 MeV) 



within a few percent. 

The correspondence rj' ■<->■ /o(980) and r] •H- /o(1500) is also characteristic of the 
"Bonn quark model" [64, 65] and the model of Nambu-Jona-Lasinio type [66] which 
include instanton interactions with axial U{1) symmetry-breaking [67]; these models 
explain the reversed mass differences between the octet and singlet states in the scalar 
and pseudoscalar nonets. 

The extra state /o(500 — 1000) is then supernumerous and is taken as glueball. This 
identification is supported at first view by the observation of this state in all "gluon rich" 
processes (i)-(v), with the possible exception of J/-?/' — )■ •jn^n^ which is difficult to select 
from background [68].^ Other results within this scheme are presented in [61] and on 
the appearance in D,Ds and B decays in [70] and [48]. The k meson is not required in 
this scheme. Further discussion follows in the subsequent sections on f^s. 

4-3. Route 3: two supernumerous states, broad glueball /o(1200 — 1600) 

A global description of processes with mesonic final states tttt, KK, rjrj, rjrf , An with 
J^*^ = 0++ in the mass range 280 -1900 MeV in terms of a X matrix and also using a 
dispersion relation method has been presented by Anisovich et al. [71] (earlier work in 
[36, 37]). Data from pp. pn annihilation into three mesons and meson pair production 
in Tip collisions have been included in a global fit. The appropriate mixing angle is 
determined for each resonance from their decays (following table 3) where the glueball 
decay is "fiavour blind" with a modification from strange quark suppression. The two 
isoscalar states in a nonet are orthogonal in fiavour space and the mixing angles of the 
respective states should fulfill 



The results of the fits suggested the classification of the isoscalar states as in table 7 with 
two qq nonets and one broad state /o(1200 — 1600) (width F = 1100 ± 140) identified 
as glueball. Alternatively, one could exchange the roles of /o(1300) and /o(1200 — 1600) 
which are both near flavour singlet. The cr meson is another extra state and it is related 
to a conflnement singularity of the qq potential. 

% Preliminary data from BES-III suggest a contribution from /o(500) centered at low masses around 
500 MeV [69]. 



m-'iQ = m2(ao) + |(m^(X*) - m^(ao)) 
2.271 = 0.970 + 1.365 = 2.335 [GeV^] 



(36) 
(37) 



^{i) _ ^ ±90°. 



(38) 
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Table 7. Spectroscopy with glueballs route 3: two supernumerous states. 

[/o(980),/o(1300)] [/o(1500),/o(1750)] nonets 

/o(12()()- KiOO) /„(.")()())/r7 glu(>1)all. <-x(ra state 



There is some similarity with route 2 in that the broad object is of gluonic 
origin. The broad "background" under the narrow resonances in elastic tttt scattering 
is represented by two poles in route 3 and by one pole in route 2 corresponding to 
the appearance of only one loop of the amplitude in the complex plane after the narrow 
states are removed. Furthermore, both schemes work without K*^{900)/k,. On the other 
hand, the fits of route 3 include the state /o(1300) which is not included in route 2. 

5. Properties of /o(500)/(7 - comparison with Kq{800)/k, 
5.1. Results on /o(500)/o' 

5.1.1. Mass and width. There are different definitions of the mass of an unstable 
particle. If P{s) ^ 1/D{s) is its propagator a conventional definition is the "on shell" 
mass with Re D{sq) = in analogy to the treatment of stable particles. A more 
fundamental definition is the "pole mass" with D{s) — 0. As this definition is a basic 
property of the (S-matrix one expects in general better properties, for example, gauge 
invariance in gauge theories; both quantities can be related (see, for example [72]). 

In the "Breit Wigner" approximation for masses around the resonance peak the 
scattering amplitude is written in terms of resonance mass tuq — ^/sQ and width F as 

^ cot 5=^^ (39) 



sq — s — imoF(s) ' moF(s) 

and the "on shell" or "Breit Wigner" mass is defined by ReT~^ = or s = sq, i-e. at 
this mass the amplitude becomes purely imaginary and the phase passes 90°; the width 
is obtained from = aX s — sq. In this way the Breit-Wigner mass and width 
are experimentally observable quantities. For /o(500 — 1000) /a one finds typically (28). 
These Breit-Wigner quantities are also relevant in the applications of the QCD sum 
rules, where they enter in the narrow width approximation to the dispersion relations as 
in (7). The complex pole mass \/i = fn — is not directly accessible from experiment 
but an analytic extrapolation of the amplitude away from the region of experimental 
measurements is required. 

A description of the observed tttt phase shifts from threshold to ~ 1400 MeV has 
been achieved using a "minimal meromorphic parameterization" [73] which is defined 
for a resonance from its complex pole position and a width F(s) oc q. Three resonances 
have been included and a background taking care of the Adler zero. Then the pole of 
the broad resonance is close to the Breit Wigner result (28). 

Another procedure to determine the pole mass is based on the application of partial 
wave dispersion relations ("Roy equations" [74]) which follow from first principles of 
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S matrix theory such as analyticity, crossing symmetry and unitarity, together with 
experimental information on nn scattering. The pole of the amplitude at complex 
energy — in — i T/2 has been found in that way as [75] 

fh„ = 441^^^ MeV, f ^ = 544^^^. (40) 

Other applications of dispersion relations for an extended range of energy {^/s < 1400 
MeV) lead to similar results with ffia = 457+i3, f ^ = 558^14 MeV [76] and rfia = 
44211, V„ = 5481}^ MeV [77] (For further discussion, see [78]). 

The large difference between the pole mass filer and the " Breit- Wigner mass" or 
"on shell mass" where the phase shift passes 90° can be understood in simple 

models such as the one by Mennessier [79] which has been slightly modified for recent 
applications [80, 81]. The analytical and unitary amplitudes for vrvr — )■ vrvr, KK are 
derived from a K matrix parameterization with a and /o(980) resonances where the 
real part of the amplitudes is derived from dispersion relations. One finds for the pole 
and "on-shell" mass 

~ 422 MeV, ~ 580 MeV; ^ ~ 900 MeV, T^^^^ ~ 1000 MeV, (41) 

so fha and are found near (40). The difference between pole and "on shell" mass 
is related in this model to the analytic properties of the denominator D{s) of the a 
propagator."*" 

5.1.2. Phenomenological appearance of /o(500)/o". In tttt elastic scattering the 
/o(500)/a contribution corresponds to a broad resonance with central mass near 1000 
MeV where the phase passes through 90°. In decays of mesons of low mass, such as 
/o(1500) — )■ (Tcr, only the low energy part of a below 1 GeV contributes because of phase 
space restrictions. In decays of B mesons with heavy b quark the broad background 
appears again which can be related to /o(500)/(7, for example in S — > Knn, B — > KKK 
(Belle [82]) and B^ — > J/t/jmr (LHCb [83]). The Belle data have been parameterized in 
[48] by a broad resonance with moving phase as in elastic scattering. There arc other 
decays, like J/ip ^ unn [84, 85] or D — )■ Stt [86] where a low mass peak around 500 MeV 
appears. The appearance of this peak in processes of relatively low energy and absence 
at higher energies may be related to a special dynamics not directly related to the a 
pole, for example to subthreshold p exchange with primary J/ip ^ pn and secondary 
p — > UJTT transition [87]. In decays like Xc ~^ crcr it is not obvious from the beginning, 
what mass spectrum the a decays will acquire. Some clue may come from a comparison 
of Xco and Xc2 with the very different phase space factors. Independently of the position 
of the mass peak, the final state tttt phase should follow the phase of elastic scattering 
if re-interactions are neglected. This was observed for J /ip ^ utttt where the S — D 
wave interference behaves smoothly over the peak region around 500 MeV as reported 
in [88]. 

+ Earlier phenomenological pole determinations by BES [85] using a variety of parameterizations gave 
the pole position = (541 ± 39) - i (252 ± 42) MeV. 
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Table 8. Ratio of couplings R„Kn = and i?^^^ = gl Jgl^ for /o(500)/a-. 



processes 






model 


reference mass 


TTTT — TTTt/KK 
TTTT — )• TTTr/KK/An 

TTTT ^ TTTT / K K / Ijlj 


0.361°:^^ 
0.561°:^, 


0.2 ± 0.04 


K matrix [89] 
S matrix [89] 

r(^K())iauc(> 1110(1(4 [90] 


a pole / ruBW 
<T pole 



5.1.3. Couplings o/ /o(500)/o' to nn, KK andrji]. The intrinsic structure of /o(500) 
can be revealed by the decay branching ratios. These can be extracted from a model 
amplitude which fits data from both channels tttt and KK. The decay ratio of couphngs 
T^aTrK — ^f'^^'^^'i^ has been derived from an extrapolation to the pole of an analytical 
amplitude or "on shell" along the real axis near tjibw ~ 1 GeV. A result of the first 
kind is obtained by Kaminski et al. [89] based on a coupled channel analysis as described 
in the previous section 5.1.1 using the analytical X-matrix model [79, 80]. They find 
faKw — 0.52 ± 0.17, i.e. a significant strange quark component in the a. In table 8 
we summarize the results on the ratios RaKn = ^"^2^^^ = \r'^r^ and Rarm which are 
derived from fits to both the elastic and inelastic tttt scattering. In an approach based 
on dispersion relations and a 3-coupled channel unitary model a similar result has been 
obtained by the same group r^i^Tr — 0.65 ± 0.18. A phenomenological description of tttt, 
KK and 7777 channels extending up to 1.6 — 2.0 GeV based on extended inelastic Breit 
Wigner formulae for /o(500)/(7 and inclusion of other resonances has been achieved by 
Bugg [90]. Near threshold the process tttt — > KK is dominated by /o(980) but a sizable 
/o(500)/(T component with RaK-K = 0.61q;2 is necessary for a good fit. We note however, 
that this quantity may be subject to additional systematic uncertainties because of 
confiicting experimental phases just above the KK threshold: the fit by Bugg follows 
the phases by Etkin et al. [91] which start rising near ipthr ~ 170° while the phases by 
Cohen et al. [92], preferred in figure 12 below, start decreasing from ipthr ~ 220°.* 

The results for R^kv in table 8 show consistency between both methods. For the 
K-matrix model it has been found [89] that TaKn evaluated for the pole mass and for 
the on-shell (BW) mass coincide. The ratios in table 8 are all derived from the same 
data and we estimate 

R^Kn = 0.5 ± 0.2. (42) 

There is also a small but significant decay into 7777 with RcrryK — 0.20 ± 0.04. 

All these results indicate that /o(500)/fT is not a pure non-strange object, but 
couples to a sizable ss (or gluonic) component. A pure fiavour singlet or glueball would 

correspond to r^jK-n = 1 or RaKn = 4/3 which is above the results in table 8. While 
some caution on these numbers is still appropriate because of possible systematic errors 
there is also the possibility of a mixing with /o(980) as discussed in [32]. According to 
these model results, a 4-quark model in its simplest form with a non-strange a [57] or 

* Bugg also obtained a consistent result R^Ktt = 0.6 ±0.1 from data on the radiative decay 
(/)(1020) — >■ 77r'^7r° with the contribution of the virtual process tttt — >■ KK within the KK loop model. 



24 



with a small mixing t^t^k < 0.1 as in [93] is not supported. 

5.1.4- Two-photon coupling o/ /o(500)/a. Following the general idea that glucballs 
couple weakly to photons one may consider the two-photon decay of the tentative 
glueball as a crucial measurement. The recent results from Belle [94, 95] on the reactions 
77 7r°7r°, tt^tt^ have renewed the interest in this problem. The two-photon width has 
been extracted by Pennington [96] as T{a — )■ 77) = (4.1 ± 0.3) keV. This result has 
been obtained from the application of partial wave dispersion relations and using the 
position of the a pole (40) from [75] . A more recent evaluation and the mean of various 
other determinations [77] is a bit lower 



A number of this order is in the range of expectations for qq states but gluonic states 
require generally smaller decay widths, for example T{a — )■ 77) ~ 0.2 — 0.3 keV [97], an 
order of magnitude smaller than what is found in (43). Therefore the glueball assignment 
has been questioned in [96]. 

A possible solution of this problem has been suggested in [80] within the analytic 
K matrix model mentioned before [79]. Given a set of coupled multi-channel 
strong processes, the amplitudes for the related electromagnetic processes are largely 
determined by unitarity and constraints from dispersion relations, although with some 
ambiguities from subtraction terms [98]. Along this general scheme, the calculation of 
the electromagnetic processes 77 — > tttt, KK, given the strong processes tttt — >■ tttt, KK 
has been carried out in [79] based on a K matrix model representing the amplitudes by 
a set of poles. The subtraction ambiguities in electromagnetic processes correspond to 
the free "direct couplings" of resonances to photons. 

At low energies the process 77 — )> tt+tt" is governed by one-pion exchange with the 
photons coupling to charged pions and this process dominates by an order of magnitude 
over 77 — > 7r°7r° where the initial photons couple to charged pions and subsequently 
re-scatter: 77 — >■ tt+tt" — >■ 7r°7r°. In the resonant scattering 77 — >■ i? — >■ tttt there is 
the possibihty of tttt re-scattering and in addition of "direct couphng" 77 — > R. At 
higher energies the re-scattering contribution decreases and the direct processes become 
dominant. They correspond to interactions where the energetic photons resolve the 
constituents inside the resonances. An example is the two-photon decay of /2(1270) 
and other tensor mesons which are well represented by the coupling of photons to the 
constituent quarks [3]. 

In [80] the "direct coupling" of /o (500) /cr has been determined from the 
simultaneous fits to the strong tttt and the 77 processes below 700 MeV. Subsequently, 
the inelastic channels including KK and a second resonance /o(980) has been included 
[89, 81]. The following results at the resonance pole for /o(500)/(7 are obtained 



r{a 77) = (2.08 ± 0.20l°:°D keV. 



(43) 



rf' ^ 77 ~ (0.13... 0.19) keV, 

^ 77 ^ (2.7 . . . 1.48) keV, V'f ^ 77 ~ (3.9 . . . 2.67) keV 



(44) 
(45) 
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where the first number comes from [80], the second one from the average of the two 
numbers in [81]. One can see that with this approach the 77 process is dominated by 
TTTT re-scattering. It is the direct couphng component which should be related to the 
constituents of the resonance and this is an order of magnitude smaller than the total 
width, quite compatible with gluonic constituents (like T ~ 0.2 — 0.3 keV, as in [97]). 
This result is based on a model consideration, but one may conclude that it is premature 
to reject the gluonic nature of a based on the 27 decay width. 

5.2. Comparing elastic tttt and Ktt scattering - a problem with the k resonance 

In the nonet classification of the scalars below 1000 MeV (route 1) the broad k and a 
particles appear. While a plays its role in all three classification routes considered the 
K appears only in route 1 and would be supernumerous otherwise if a real resonance. In 
this subsection we add a few comments on these broad objects in a comparison. 

At first sight there is a similarity in low energy tttt and K71 scattering and the 
strong low energy interaction suggests the existence of two broad resonances a and 
K with a width comparable to mass. The first narrow resonance of low mass in tttt 
scattering appears with /o(980) and in Kn scattering with Kq{1A30). A closer look 
at the scattering phase shifts, however, reveals a considerable difference between both 
scattering channels. 

In the elastic scattering region we may parameterize the amplitude for the 
production of a resonance R above a broad background B, following Michael [99], as 
product of the respective S matrices which satisfies unitarity (also called "Dahtz-Tuan 
method" ) 

S^SrSb, or T = TB + e2'^«TR (46) 

where S — 1 + 2iT and T — sin Se^^. This approximation applies below the inelastic 
thresholds of tttt KK at m^^ ~ 1000 McV and Kn Kt]' at tjik^ ^ 1450 MeV 
respectively (disregarding the negligible Krj production) and then the overall phase of 
S = e^^^ in 5 = 5r + 5b- According to [99] (46) is also applicable for inelastic resonances 
and near elastic backgrounds. If the contribution of the lowest narrow resonances R to 
the phase 5 are "removed" the remaining background phase 5b shows quite a different 
behaviour for the two channels (sec figure 5): 

m^^ ^ 1000 MeV : 5b = 90° (47) 
niK^ ^ 1400 MeV : 5b < 40°, (48) 

i.e. in tttt scattering the "background" itself shows a resonant behaviour whereas in Kn 
scattering the background phase saturates already below 40° according to the present 
experimental status of Ktt scattering which presents itself as follows: 

There is still a discussion on the appropriate description of low energy Ktt 
scattering. Elastic Ktt scattering amplitudes have been reconstructed from processes 
K^p — )■ K^TT'^n and K^p — >■ X^tt^A"'""'" by isolating the respective One-Pion- Exchange 
contribution [106, 105]. Following the LASS Collaboration [105] the S-wave contribution 
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Figure 5. Phase shifts of elastic scattering, Left: tttt phase shifts from Dispersion 
Relations [100] and from resonance fit to CM-I data [101] (/o(980) + background 
[102]); Middle: Ktt phase shifts (i4:o*(1430) + background [103, 104]); Right: Ktt phase 
shifts from BABAR [103] {D+ R-n+e+v^), LASS [105] (if^p ^ K^^ir+n) and 
Estabrooks et al. [106] {K^p K^tt-A++), figure from [103]. The Ktt background 
phase stays below 40°. 



with I = \ can be represented by the superposition as in (46) of the i^'o(1430) Breit- 
Wigner amphtude and a background which is parameterized by an effective range 
formula (cms momentum k) 

, 1 fop 
cot(5BG) = ^ + ^; (49) 

a recent reanalysis has been presented by the BABAR Collaboration [103, 104] which 
gave the numbers a = 1.95 ± 0.09 GeV^^ and b = 1.76 ± 0.36 GeV^^ for the / = | 
S-wave. 

A model independent determination of the Kn phases in the elastic region is 
possible from the decays of D — )■ KneUe and r — )■ Kt^v^- as in these decays, according to 
the Watson theorem [107], the Ktt production phase agrees (modulo tt) with the elastic 
Kn scattering phases. 

A measurement of all the decay distributions in D — )■ Kireu^ by the BABAR 
Collaboration [103] has determined the I = \ Kn S-wave phase shifts in the range 
700 — 1400 MeV. The results agree nicely with the earlier results on Ktt phases from the 
production off the proton target and confirm the description of the S-wave amplitude 
in terms of the i^^^(1430) and the effective range background (49) (see figure 5, right 
panel). The resulting background phase approaches 38° at around 1400 MeV where it 
reaches a maximum, jj 

An analysis of r — KnVj. decays has been presented by Belle [109]. The Kti mass 
spectrum requires contributions in addition to i^*(892) and a successful fit includes an 
S'-wave with a superposition of Kq{^QQ)/k and i^Q(1430). Also a LASS-type fit with 
background (49) was possible but found only marginally consistent with previous values 

jj The existence of an S'-wave with a phase variation consistent with (49) has been noted before [108] 
in a limited mass range around i^*(892) with lower statistics. 
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for the parameters a and b. However, the measurement of further decay observables 
is required for an independent determination of the Ktt phase shifts and a test of the 
Watson theorem, to compare with the resuhs from D — > KneUe- First results have 
been presented by BABAR with similar fits [110]. At present we stay with the results 
presented in figure 5. 

The situation with a low mass resonance Kq{800)/k, is discussed by the PDG [78]. 
The existence of this state close to the Kn threshold and with a large width around 
500 MeV is considered difficult to estabhsh. A positive result is obtained recently in an 
approach based on a dispersive representation of the Ktt scattering amplitude [111]. 
After continuation into the complex energy plane a pole has been identified at 

m« = 658 ± 13 MeV, r« = 557 ± 24 MeV, (50) 

just above Kn threshold at around 640 MeV. 

In summary, descriptions of the scalar Kn amplitude with a k pole built in are 
possible but the apparent slow phase motion of the corresponding amplitude which 
saturates below 40° is untypical for a resonance. In this way a characteristic condition 
for the existence of a resonance - the phase shifts arc passing through 90° - is not fulfilled. 
One could think of a repulsive background under the k, with destructive interference 
[112], but then the background phase ought to be larger than the resonance phase 
{\Sb\ > \^r\)- Alternatively, an enhanced but not necessarily resonant low energy Ktt 
interaction could reflect the strong t-channel tttt — > KK low energy (a) forces. ft 

6. Intrinsic structure of /o(980) 

6.1. Mixing angle from various observables 

We discuss here several observables within the qq constituent model and determine 
the qq mixing angle defined as in (31) following earlier studies [61, 70]. Some of these 
observables have been also considered in [113]; other determinations of the mixing angle, 
also in the context of a 4-quark picture, have been reported in [114], but no clear picture 
has emerged. We consider the following ratios of branching ratios: 

(i) Decay into tttt and KK. This decay is often described in terms of a Breit Wigner 
resonance formula, following Flatte [115], with amplitude 

y ^ 9i9j^/PIPj 

so-s- i{pngl + Pk9k) 

where gi denotes the couplings of the resonance to the tttt and KK channels and pi = ^ 
the respective phase space with the momenta ki in the resonance rest frame. The 
BES Collaboration has measured the decays J/ip 0(lO2O)7r^7r^/0(lO2O)i^"'"i^~ and 
determined gj^/gl = 4.21 ± 0.33 (errors added in quadrature) [116]. 

ft This viewpoint has been envisaged in discussions with P. Minkowski. 



(51) 
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Various results have been obtained from data on the process tttt — )■ KK together 
with TTTT TTTT in an apphcation of different coupled channel approaches. The data 
mainly used in these analyses because of their accuracy are from Cohen et al. on 
K^K~ [92] and from Etkin et al. [91] on K^K^ final states. In table 9 we list some 
results on glc/gl which included KK data in the fits: from Martin et al. [117] using 



Table 9. Couplings gj^/g^ in (51) for /o(980) from tttt — > KK and from J/tp decay. 

Martin Estabrooks Bugg Kaminski Mennessier mean (tttt) BES {J/i>) 
4.0 ±0.6 3.2 ± 1.8 2.3 3.9 ±0.9 9t" 3.3 ± 0.6 4.21 ± 0.33 



several different models for resonances; from Estabrooks [118] based on a K matrix fit; 
from Bugg et al. [119] using a global parameterization based on a K matrix, see also 
[113] for interpretation; from Kaminski et al. [120] using a unitary model with separable 
interactions (see also [89]) and the analytic K-matrix model by Mennessier et al. [81]. 
As the ratio Qk/qI depends essentially on the same KK data we average the weighted 
results as usual but assume a full correlation of the input data and obtain the mean in 
table 9. Finally, together with the independent BES measurement we find 

/o(980) : gl/gl = 4.0 ± 0.3; gx/g. > 0. (52) 

The same sign of the couphngs gx and g^^ has been derived from the fits [117, 118] 
to the K~^K~ data by Cohen et al. In these data the relative sign of the 7r~^7i~ — )■ K^K~ 
and Ti^Ti~ — )■ TT^TT~ amplitudes has been fixed from the observation that near 1 GeV the 
P wave in both channels is governed by the tail of the p and assuming these amplitudes 
to be in phase. This observation resolves the previously noted sign ambiguity (for 
example in [113, 81]). The result (52) can be expressed in terms of the mixing angle (psc 
defined in (31). Using table 3 with (f)sc = a + 7r/2 according to the definition (15) we 
find g\/gl = ll^Ml^ ^ 

gl/gl^{l + V2cotct>scf/2> (53) 

in the symmetry limit p = 1 and yields <j)sc — 29.9° ± 1.3° (for gK/g-K > 0, otherwise 
(j)sc = 162.4° ± 0.5° noting that gx+K- / g^n+^K- = (1 + v^cot0sc)/2 from table 3). 

(a) Decays J/%1) ^ /o (980)0 and /o(980)a;. In these decays the ss and nn quark 
components of /o(980) are projected out and according to (31) their ratio can be related 
to the quark mixing angle [61] assuming the dominance of the "singly disconnected 
diagram" [40] 

^'/^ ~ B{J/i^ ^ /o(980)a;) " ^'1^ ~ 

with cms decay momenta ki of (f) and uo. Taking these branching ratios from PDG as 
in (30) we find cot^ 0,^ = 2.46 ± 1.12 or = ±(33 ± 7)°. This result is similar to 
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the corresponding ratio for 77', -Rj/^ = 2.05 ± 0.43 using PDG values and this yields 
4>ps{il') = (35 ± 3)° to be compared with the experimental mixing angle (j)ps = 42° in 
(17a). This observation was the original motivation for the mixing ansatz for scalar 
mesons in [61]. 

(Hi) Two photon decays o//o(980) and ao(980). These decays are sensitive to the 
charged constituents of these mesons. In an update of [61] we consider the ratio R^j in 
a qq constituent model 

r^^(fo(980))' r^^^'E^^^?)' (55) 

with amplitudes Cq for q = {u, d, s), (X^c^ = 1), for the contributions from qq of charge 
Qq and obtain with mixing (31) 

i?/o^o ^ x^cos0,e)'. (56) 

The two photon decays have been studied recently by the Belle experiment for /o(980) 
[94, 95] and ao(980) [121]. We take here T^o^ = 0.29^[j:|^^ keV for /o(980) from the 
PDG fit, in case of we use their branching ratio F^'^/F^^ = 0.183 ± 0.024 to obtain 
5(ao(980) ^ Tvri) = 0.85 ± 0.02 and then F^o = 0.25lg:[j^ keV; finally, with (56) 

= 1.16 ±0.41, 0,ei = (23.7 ±8.1)°; 0,e2 = (125.8 ± 8.1)°. (57) 

(iv) Charmed meson decays D,Ds — >■ /o(980)7r. The decays of D'^ and Df into 
/o7r+ and Xo(1430)7r+ provide a further possibihty to determine the flavour structure 
of /o(980) and we follow here the investigation in [70]. They consider the two ratios 

_ r{Dt ^ /o7r+) _ r(Dt ^ /o7r+) 

™ r(L>+ ^ /o7r+) ' r{D+ ^ K*{U30)7r+) ^ ' 

Assuming again the flavour composition (31) for /o(980) and the iro(1430) to belong to 
the same nonet, we write the amplitudes for the three processes involved as follows 

(a) Dt ^ /o7r+ A = cos 0,cKdK> (59) 

sin 1 



(6) ^ /o7r+ B 



-VudV:Al + e) + COS0,eK.K> 



a (60) 



(c) ^ i^*(1430)7r+ C = KdV;:(l + e)a. (61) 

Process (a) is given by the amplitude with Cabibbo- favoured decay c — >■ sW^ , 
— )> M(i — )> 7r+ followed by the spectator interaction ss — )■ /o with strength a. In 
the Cabibbo-supprcsscd process (b) there is the same direct 7r-emission diagram as in 
(a) but one has to add the "colour suppressed" amplitudes with decay c ^ d W^, 
W'^ — >■ ud, dd — >■ /o and c — )■ s W'^, W'^ — >■ us, ss — >■ /o both of relative size e. Finally 
in (c) there is one direct and one colour suppressed amplitude. In the approximation 
Vud = Vcs = cosi?c, Vus = -Vcd = sin-dc, with Cabibbo angle -de 

Rdi = '2- coe cot^ 0,,- ^ — - (62) 

?6 |1 - (V2cot0sc - l)e|^ 

Rd2 = — COS^ <t>scTT—^ (63) 
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where qi is the phase space for S wave decays with pion momenta in the rest frame of 
processes (a — c). Using the branching fractions estabhshed by the E791 Collaboration 
[86] the following results have been found [70] 

cot2 05c/|l - (v^cot(/),c- l)ep = 1.26 (1.0 ±0.4) (64) 
cos^(/),c/|l + e|^ = 0.52 (1.0 ±0.3) (65) 

with the two solutions 

cot <Psci = l.lllai, = (42.0ty)°, ei = (2.85 ± 5.35) IQ-^ (66) 

tan0,,2 = -0.34^0;^6, = 161.2;?|D°, ^2 = 0.31;°;°°^ (67) 

It is satisfactory that the colour suppressed amplitudes are indeed found small. 

Table 10. Quark structure of /o(980) from four observables, mixing angle as in (31). 



observable exp. result sol. (psci sol. (psc2 ry'-like as in (35) 



91/91 


4.0 ±0.3 


(29.9 ± 1.3)° 


[(162.4 ±0.5)°] 


3 


T3<t>/'^ 


gK/g^T > 


0° < (psc < 90° 






2.29 ± 1.04 


(33 ± 7)° 


(147 ±7)° 


2 




1.16 ±0.41 


(23.7 ±8.1)° 


(125.8 ±8.1)° 


f = 1.815 


Rm , Rd2 


(64),(65) 


(42.0lf j)° 


i6i.2i?i:D° 





6.2. Final result on mixing angle 

Our results are summarized in table 10. There is a good agreement between the four 
determinations for mixing angle solutions with (psci < 90°, contrary to (I)sc2 > 90°. 
This can be taken as an independent confirmation of the relative sign of gx/ g-n in (52). 
Finally, we obtain as average = (30.2 ± 1.2)° which is dominated by the observable 
glklgi- In view of systematic uncertainties, such as strange quark suppression effects, 
we estimate the total error to be around 10%, therefore finally 

/o(980) : <Psc = (30 ± 3)°. (68) 

The mixing angle for /o(980) is then near the mixing angle (17&) for 77', i.e. (psc ~ 
(pps = 42° as in route 2, or approximately with = arccos ^273 = 35.3° for 
|/o(980)) (1,1,2)/V6 as in (35) (observables in table 10 ''r]'-\ike") and so /o(980) 
is close to a flavour singlet state; octet would have (f)^^ > 90°. One of the two solutions 
found in [113], 23° < 4>sc < 37°, is consistent with our result (68). 

6.3. Possible gluonic contribution 

While the data are nicely consistent with /o(980) being a pure quark state, we may 
estimate a possible gluonic contribution as follows. The fit result (68) which is largely 
dominated by the result on gj^/gl corresponds to a value R^j = 1.53 according to (56), 
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so the actually measured value in (57) is about 75% of the nominal value. Taking into 
account a gluonic component in an extension of (31) by 

|/o(980)) — cos l3 (sin^scl^^^^) +cos0sc|ss)) + sin /3\gg), (69) 

we can write a properly extended equation for T^^ in (55) where the additional gluonic 
term is of 0{al). Neglecting this term here we obtain with the above result 

R^^ ^ Rf^ cos^ P; cos^ p ^ 0.75 ± 0.25, 

i.e. the fraction of a possible gluonic component is about ^ ± |. Note that the addition 
of a gluonic component would leave R^"!'^^ unchanged for fixed qq fractions but would 
modify the other two observablcs considered. 

Although with a considerable uncertainty, this result suggests that more precise 
data on the above observables could further strengthen the picture and provide better 
estimates on the amount of glue inside this meson. Recent measurements of decays 
Dg — > Sn and Dg — > KKn [122, 123, 124] have shown the potential of independent 
determinations of particular, BaBar has compared the Tin and KK mass 

spectra just around 1 GeV and found good agreement in their shape [123] up to about 
1.1 GeV which should directly determine this observable. 

More precise results on Rjf^^^^ and ri^^j^^^^ could improve our knowledge of the 
flavour structure of /o(980) and arc accessible by ongoing or new experiments. A word 
of caution on the interpretation of the 77 results is in order. The model calculations 
suggested a considerable contribution from hadronic re-scattering in the cr ^ 77 decay 
(section 5.1.4). On the other hand, the corresponding calculations for /2(1270) have 
shown the dominance of the "direct" contribution as expected [81]. Results on /o(980) 
in between the extremes indicate re-scattering contributions but no strong difference 
between direct and PDG values for /o(980) — > 77; a strong model dependence on 
the contributing exchanges can be observed [81]. Further studies, for example with 
dispersive representations, could strengthen the conclusions. In this presentation we 
assume the dominance of the direct contribution in (56) as for /2(1270). 

7. Experimental evidence for /o(1370) and /o(1500) - a reassessment 

A crucial precondition of the mixing scenario in route 1 with a supernumcroTis state in 
the mass region around 1600 MeV is the existence of /o(1370). This section is devoted 
to a critical reassessment of the evidence for this state in various processes and the 
comparison with /o(1500) nearby in mass. In previous surveys supportive [125] and 
sceptical views [5] have been presented. 

7.1. Production in pp annihilation 

The state /o(1370) has been discovered in pp annihilation at rest by the Crystal Barrel 
Collaboration (CBAR) and has been studied together with /o(1500) in the reactions 
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Figure 6. Upper panel: Dalitz plots for pp — ^ Stt^ annihilation at rest; Bottom: 
Dalitz plot for pp — > n'^r]!] together with Mass^ (7^77) spectrum from Crystal Barrel 
Collaboration [126] (figures from [5]). 



[126, 127, 49] 

a) pp — i- Tr'^Tr'^Tr'^, b) pp — i- 7r°?7?7, c) pp — > TT^vr'^r]. (70) 

The analysis is based on the Dahtz plots shown in figure 6. The Svr" Dalitz plot has 
sixfold symmetry and is based on ~ 700000 events. The coloured version of the top 
plot shows events in bands near the corners at m^(7r°7r°) ~ 2.25 (= 1.5^) GeV^ which 
corresponds to the isotropic decay of the S wave resonance /o(1500) (at fixed (vertical) 
mass squared S23 the (horizontal) S12 depends linearly on the decay angle cos 6*23). At 
the mass near = 1.370^ = 1.9 GeV^ there is no band of comparable strength, the 
two bumps at the edges of the plot correspond to /2(1270) with its peaked angular 
distribution according to spin 2. The evidence for the second resonance /o(1370) is 
based on the results of a multi-resonance K matrix fit. 

A more direct evidence is suggested from reaction (70b). In the n^rirj Dalitz plot 
(Bottom left of figure 6) one observes the horizontal and vertical bands from ao(980), in 
the diagonal along the arrow "glueball" there is again a band with fiat density related 
to /o(1500) — )■ rjrj. In parallel, a second band appears to the right (at lower rir] mass). 
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Figure 7. Upper part left: Dalitz plot for pp — )• Tr^r^ry at ^/s = 2.0 GeV (Crystal 
Barrel); right: mass^ {'qr/) spectrum (figures from [128]). Bottom: 7777 mass spectrum 
in pp — >■ -K^-qr] at ^/s — 3.0 GeV, figure from [129] (Experiment E-760). 

but the band does not continue to the edges of the plot as the /o(1500) band does. 
The projection in m^(?7?7) (Bottom right of figure 6) shows peaks at higher mass with 
contributions from /o(1500) and at the lower mass near m ~ a/1-7 = 1.3 GeV, and this 
peak has been related to /o(1370). The overlap with the ao(980) bands makes it difficult 
to directly judge the relevance of these peaks and one relies on the fitting procedure. 

A clarification can be expected from a change of primary energy which would 
shift the diagonal bands at fixed r^rj mass away from the crossing r/vr^ bands. Such a 
measurement has been performed with 900 MeV antiprotons {^/s^ = 2.05 GeV) again 
by CBAR experiment [128]. The rjrjir^ Dalitz plot with two- fold symmetry is shown in 
figure 7. Indeed, the diagonal band (C) which corresponds to /o(1500) is now shifted 
away from the crossing ao(980) bands (A). The new bands (B) are related to a2(1320). 
A second diagonal band, between /o(1500) and the crossing point of the 02 (1320) bands 
would be expected for /o(1370), but there is no indication of any second band. Also 
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the peak in the rji] mass spectrum has changed into a shoulder easily related to the 
resonances in the rjn^ channel. It is difficult to invent a mechanism which would make 
disappear a strong /o(1370) signal visible at ^/s^ = ^rUp by increasing the cms energy 
by only 0.17 GeV. 

At yet higher energy the process pp tt^tji] has been measured by Experiment 
E760 at Fermilab [129]. At ^/s^ ~ 3 GeV the mass spectrum in figure 7 (lower part) 
again shows a clear peak at 1500 MeV suggestive of /o(1500) (see also discussion in [5]), 
but the shoulder at low mass has now disappeared entirely and there is no indication of 
/o(1370). More recently, Experiment E835 has analyzed data on the same reaction at 
^/s^ = 3.4 GeV and determined the spins of the observed resonances [130]. As a result 
they identified their lowest mass peak with /o(1500) and their second peak as /o(1710). 

At all energies in these reactions /o(1500) shows up clearly as band in the Dalitz 
plot and sometimes as peak in mass spectra. On the other hand, no such "direct" 
evidence is visible for /o(1370). The effect appears only as a result of a global resonance 
fit at the lowest primary energy where there is a large overlap with other resonance 
phenomena. 

7.2. Search in phase shift analysis o/tttt scattering 

A method which does not rely on parametric global fits is the determination of the 
complex scattering amplitude in steps of increasing energy directly from the data. Any 
resonance decaying into the tttt final state has to appear also in the elastic tttt scattering. 
The scattering amphtude T^is) — xmor{s)/{so — s — imor(s)) with a Breit-Wigner 
resonance pole at sq — irrioT moves along a circle for increasing (real) mass m = \/s 
around (0, ix/2) with radius ,t/2 in the complex plane ( "Argand diagram" ). This circular 
behaviour is smoothly modified by background effects and it is a characteristic signature 
of resonant behaviour to look for. The aim of phase shift analysis is the determination 
of the complex scattering amplitude from the scattering data. 




n,a, .O2 



p N,i 

Figure 8. One-Pion-Exchange: tttt scattering contributes to irp — )• -k-kN{^) reaction. 



7.2.1. Methods to reconstruct tttt scattering amplitudes. In this subsection we will recall 
the main assumptions underlying the reconstruction of the tttt partial wave amplitudes 
from measurements of the reaction 

7r~p — >■ 7r~7r''"n (71) 
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(and others) using the One-Pion- Exchange model (OPE) [131] (see figure 8). The 
incoming pion scatters off a pion from the pion cloud surrounding the target proton 
at small momentum transfers t = {p^^ — Pout)'^- The theoretical description of the 
production process is based on Regge theory and the first aim is to isolate the one- 
pion exchange contribution from other exchange contributions (for a review, see [132]). 
The hehcity amplitudes for process (71) in the "t-channel" (Gottfried Jackson) frame 
(quantization axis along exchanged pion) obtain contributions 



fix, = N^^^^ Te 5,,^y {m^.^^Ao + ^ VW+i) <^n,o) + ■ ■ ■ (72) 

with the TTTT — >■ TTTT partial wave amplitude = sin^^e*''* as a factor; A'" is a known 
normalization, the TTTT mass and q = — . The first term with 

the pole at t = corresponds to the OPE Born-term with the pion pair produced 
with helicity fi = 0. The nucleon helicity flip amplitudes with A = —A' dominate 
over non-flip amplitudes at high energies and vanish like ^Z^. Simple models suggest 
corrections from absorption at small impact parameters of the collision for net helicity 
flip n — \X — X' + n\ = 0, such as the second term in (72) [133, 134, 135, 136, 137]. 
Furthermore, contributions from 02 exchange are found important at larger t while Oi- 
exchange had been neglected at the time. The determination of the tttt amphtudes 7) 
can be based on (72) (see, for example [136]). 

Alternatively, one may use t-integrated angular distributions W {cos 6 , (p) , 
represented by the spherical harmonic moments {Y^'^) as function of uIt^t: which are 
expressed by the density matrix p and helicity amplitudes [137] 



(1?) = E '^Z.'^pZ.' (74) 



AA' 

It is useful to introduce amplitudes and the density matrix for asymptotic natural (-I-) 
and unnatural exchange (— ) which do not interfere in the observables considered here 
using the parity relations fl^,^ = (-)^+^'+''/1a-a'-m 

9SI -lifxx'.Ti-rfiy-,) (76) 

fC^' - \{C' T {-ff^:^) (77) 

In general, there are more density matrix elements than measured moments. Motivated 
by the models for OPE with absorption and amplitudes as in (72) some simplified 
assumptions for the structure of density matrix elements have been proposed [133]: 

(i) "Spin coherence:" There are only nucleon flip amplitudes. In case of experiments 
with unpolarized target the same consequences on observables follow already from the 
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proportionality between flip and non-flip amplitudes 

9r.l = -'^Vn. (78) 

This relation also implies for the density matrices Rank (p*-^^) < 1 while Rank (p*^^^) < 2 
would be allowed for in general. Then one can also write 

with gf!:'^'' = y^l + At this point the sum over nuclcon hclicitics has 

disappeared and the relation between the observable moments and the nn production 
amphtudes is very similar to on shell tttt scattering and one uses for the amphtudes g^^ 
the spectroscopic notation L^,^^ — Sq~\ Po~\ Pi^^ ■ ■ ■ 

(ii) "Phase coherence:" A second simpliflcation suggested by absorbed OPE (72) is 
the appearance of the same phase 5^ for all helicity amplitudes with the same tttt spin i. 

These two assumptions are motivated only for small t; they could be violated by 
the presence of other exchanges like ai or a2 exchanges. Both assumptions together 
actually lead to constraints between the observable moments (V/^) which arc rather 
well satisfied experimentally at 17 GeV up collisions [101]. Also the reduction in rank 
of the density matrix has been confirmed experimentally near the p meson mass peak 
[138]. 

Alternatively, one can relax the assumption on phase coherence and introduce extra 
phases between the different helicity amplitudes [134] . Measurements of relative phases, 
such as = (^1^ ^ — (^0^ ^ , can check the quality of the phase coherence approximation 
and the dominance of a single production mechanism for spin £. Further discussions of 
the two-pseudoscalar system and the ambiguity problems are given in [139]. 

As an example we give the relations for identical particles with even angular 
momentum [101, 139]. 

^{yi) = "fmrn cos(0^- - 0^-) + v^i^o^ipri cos(05- - 0^-) 

MY,')^l\D^\'-lm\' + \Dt\^) 

Va^{Y!) = ^[L'o ll^rl cos(0^- - 0^-) (80) 

We stress again, that this kind of equations without reference to nucleon helicities rely 
on the assumption of "spin coherence" for unpolarized target experiments. 

The vahdity of the results obtained from production experiments based on the OPE 
model has been demonstrated recently by the nice agreement between Ktt phase shifts 
obtained from Kp — 7> Kim and from decays D — >■ KireVe using the Watson theorem (see 
subsection 5.2 and figure 5). 

7.2.2. Results on 'k^tt~ — )■ 7r"'"7r~. A measurement of Ti~p — )■ TT~Ti^n with unpolarized 
target has been carried out by the CERN-Munich Collaboration [137]. They measured 
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the moments (Y^) in the range 0.6 < 771,^71 < 1.8 GeV from which the partial waves 
with TTTT spin £ < 3 have been determined. The overall phase of the amplitude cannot 
be determined directly but it is inferred from the Breit-Wigner phase of the leading 
resonances p(770), /2(1270) and p3(1690). In general, the moments with L < L^ax 
determine the partial wave amplitudes only up to discrete ambiguities. Writing the 
amplitude T{s^z) = f{s)Il-I^{^''' {z — Zi) as function oi z — cosO one notes that 
the measurement of \T{s,z)\'^ at fixed s — rri^^ does not determine the sign of the 
imaginary part of the zeros Zi ("Barrelet zero" [140]) and so one obtains 2^"""^/^ different 
solutions. Some can be removed by unitarity constraints. Subsequently, an experiment 
with polarized target has been carried out by the CERN-Krakow-Munich Collaboration 
[141, 142]. The following results have been presented on tttt scattering amplitudes 
above 1 GeV which are of interest for our discussion of /o(1370). 

(i) The first analysis by Hyams et al. [101] determined the moments averaged over 
\t\ < 0.15 GeV (CM-I). These moments have been fitted by amplitudes 5'^''^^ (/x = 0, 1) 
as in (73)-(79) parameterizing the tttt amplitudes as function of m^^r by a if matrix. In 
addition, local deviations from this solution have been studied. The S wave featured 
the narrow /o(980) interfering with a broad background passing 90° near tUt^t^ ^ 1 GeV. 
A classification in terms of amplitude zeros has been given but no study of alternative 
solutions has been attempted. 

(ii) The first full study of ambiguities from CM-1 data has been carried out by 
Estabrooks and Martin [143, 144]. They obtained four different phase shift solutions 
A-D classified according to the Barrelet zeros. A K matrix fit to the S waves of the 
four solutions and taking into account KK data as well, yielded three scalar resonances 
/o(800), /o(1005) and /o(1540) with width around 1000, 8 and 200 MeV resp. [118]. 
A subsequent analysis in which the phase of the amplitude has been determined using 
the analyticity constraints of dispersion relations produced modified solutions a, ^ {(3') 
starting from A,B [145]. Because of earher 7r°7r° data the solution C as starting solution 
has been rejected at the time but has been favoured later. 

(iii) A second analysis by the CERN-Munich group [146] (CM-11) included in the 
fits to the primary data a complete error correlation matrix using larger mass bins and 
obtained in that way tttt amplitudes more strongly constrained than those of [101]. Also 
in these fits some t dependence is parameterized. The ambiguities are analyzed and 
again 4 solutions are identified. The smaller error bars of zeros Im Zi do not allow the 
crossing of Im zi at 1.5 GeV as in the previous analysis (i) in [101]; therefore these 
earher results should not be trusted above ~ 1.4 GeV. 

(iv) A further experiment with transversely polarized target has been performed 
by the GERN-Krakow- Munich GoUaboration [141, 142] with the aim to check the 
assumptions underlying the previous experiments. Large polarization effects have been 
found which prove the existence of sizable nucleon non-ffip amplitudes of ai-exchange 
type, so far neglected. However, the more general property of "spin coherence" is still 
approximately valid, i.e. amplitudes with nucleon spin flip and non-flip are proportional 
for different di-meson spin orientations £, /ito a. good approximation in the average over 
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Figure 9. Dependence of S wave intensity on 7r°7r°-mass in Tr~p tt^tt^tt,. The first 
two dips reflect the destructive interference of a broad background with the narrow 
/o(980) and /o(1500), figure from [150] (GAMS Collaboration). 

t as in (78), so that the previous analyses can be maintained (see also discussion in 
[61]). In [147] it has been shown that the tttt S wave phases in the region 1000-1400 
MeV agree quite well from the polarized target experiment [148], from analysis (i) [101] 

and solutions ( ) or ( — I — ) from [146], corresponding to solutions A,C in [144]. 

We therefore consider only these two solutions in the following. 

7.2.3. Results on tt+tt" — tt^tv^. Partial wave amplitudes of the charge exchange 
reaction have been obtained from the production process n^p n^n^n at incident 
momentum 18.3 GeV by experiment E853 at the Brookhaven National Lab [149] and at 
100 GeV by the GAMS Collaboration at CERN [150]. Both experiments measured the 
moments integrated over small \t\ and determined the partial wave amplitudes using 
methods as discussed above. The E853 experiment found for mTTTr 1-5 GeV one 
solution with phase coherence {(t>Do — 4>D- ~ 0) and another one with large relative 
phase; above that mass the errors are rather large. Also at larger —t > 0.4 GeV^ the 
phase coherent solutions disappeared. GAMS found solutions with phase coherence in 
the full mass range studied up to 3 GeV and imposed this condition subsequently in the 
phase shift analysis. The mass spectrum of the S wave is shown in figure 9. 

7.2.4- Selection of a unique isoscalar tttt S wave. A further selection between the two 

phase shift solutions ( ) and ( — I — ) or A and G, respectively, left over from 

the TT'^TT" analysis can be obtained by comparing to tt'^tt" — >■ tt^^tt^ with its different 
ambiguity structure. Comparing GERN-Munich and GAMS data which have good 
accuracy up to 171-,^^^ < 1.8 GeV a unique solution has been estabhshed [102, 151]. J A 

I These analyses are preliminary concerning the parametric S matrix fits. 
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Figure 10. Argand diagrams for tttt Sj waves: Sq solutions ( 1 — ) and ( ) 

from TT+TT" — 7r+7r~ (CM-II) but using new inelastic ^2 as shown; right panel: 6*00 
amplitude for tt+tt" 7r°7r° (GAMS) and component thereof (figures from [102]). 

new (inelastic) 1 — 2 amplitude S2 for the mass range 1000 < 171^^^^ < 1800 MeV has 
been constructed from data [152, 153] and has been used to obtain updated 7 = 
amplitudes So (see Argand diagrams left in figure 10 for new Sq and S2). The diagram 
rightmost in figure 10 shows the S wave amplitude ^oo for the process 7r"'"7r~ — >■ tt^^tt^ 
which has been reconstructed from the magnitude |<5'ooP, in the given normalization, 
the phase \(j)s — 4>Do\ &nd the Breit Wigner phase and normalization of /2(1270) in the 
-Do wave. The / = amplitude So is obtained from 5*00 = So — S2 and is shown in the 
same diagram (normalization \Sj — || < |). 

The amplitude Sq obtained from tt^tt" data should agree with 5*0 from n~^7T~ . This 
is only true approximately for the ( — + — ) solution: there is a circular motion from 1000 

to 1740 MeV with a smaller circle superimposed above 1400 MeV. Solution ( ) 

would turn into the center of the circle instead and so it has been rejected. There are 
discrepancies between the remaining solutions Sq which have been related to systematic 
uncertainties in the overall phase, the 1 — 2 amplitude and influence from higher partial 
waves. These uncertainties over large mass scales should not affect the nature of the 
local resonance phenomena like /o(1370) and /o(1500). 

The amount of agreement is nevertheless impressive, in particular concerning the 
existence of /o(1500) which shows up clearly as a small circle in both analyses. Mass and 
width are roughly compatible with the PDG values [4] (m=1505 MeV, r=109 MeV). 
The partial width x^^r — ^mr/^tat in the resonant amplitude Tq — XT^T^moT/ {so—s—imoT) 
can be estimated from the depth of the /o(1500) circle in figure 10 (CM-II), see also 
inelasticity 77° in figure 11 (Ax = 2A770): 

/o(1500) : x^^ = 0.25 ± 0.05 (CM-II), x^^ = 0.349 ± 0.023 (PDG). (81) 

Here the first result (CM-II) is determined from Im Tq of the resonant elastic partial 
wave amplitude and the second one (PDG) from all inelastic channel cross sections; 
both should agree because of the optical theorem and they roughly do within 30 %. 

In order to represent the selected solution So{ — h— ) in terms of a minimal number of 
resonances an ansatz for a unitary 3-resonance S matrix for the 3 channels (tttt, KK, An) 
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Figure 11. Resonance fit (83) in comparison with data (CM-I/II): Argand diagram 
for corrected So wave, phase shift 5q and inelasticity 7]^ (figures from [102]). 

has been fitted to the CM-II data. In a generahzation of (46) the S matrix is written 
as product of three unitary matrices for resonances Sr — 1 + 2iTR 

^ ~ 'S'/o(980)'5'/o(1500)'S'broad (82) 

Tr = [ml - ml^ - i{pigl + p2gl + Pzol)]'^ 

X P^'^{9i93)p^ (83) 
where pi — Ikij \fs and {gigj) is the matrix with real constant couphngs. This S 
matrix behaves locally like a superposition of a narrow inelastic resonance over a largely 
elastic background as in (46), but it docs not fulfill the symmetry requirements for the 
S matrix exactly, therefore further tests are necessary. Also, the formulae (82), (83) 
do not represent the correct behaviour near tttt threshold and cannot be used for the 
determination of the pole of S'broad- 

These formulae have been fitted in [102] to the data from CM-II for mT^-^ > 1000 
MeV and the phases from CM-I in 600 < m^7r < 1000 MeV for definiteness. As can be 
seen in figure 11 the three resonances give a good overall description of the data. The 
broad state corresponds to a Breit-Wigner resonance with mass parameter mo = 1100 
MeV and a total width of similar size which we relate to the high mass tail of /o(500). 

7.2.5. Search for /o(1370) in tttt scattering. The fit in figure 11 reproduces the main 
features related to the rapid phase variations for the two narrow resonances at 1000 and 
1500 MeV above a broad background which can also be described as resonant state. In 
particular, there is no indication of an additional state around 1370 MeV which should 
become noticeable as an extra circle and the related variations in phase and inelasticity 
comparable to /o(1500). If we estimate from figure 11 that an extra state with more 
than 1/3 of elasticity x^tt of /o(1500) should have been noticed, we arrive at the limit 

/o(1370) : < 0.1 (CL = 95%) (CM-II data) (84) 

The Argand diagram of solution C in [143] shows a similar behaviour. The K matrix 
fits to those data require /o(1500) but no additional resonance near 1300 MeV [118]. 
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It should be noted that in figure 11 the data points represent in each mass interval, 
together with the results from the other partial waves an (almost) perfect fit to the 
original (Yl^) data at this mass. There is no correlation between data points at 
different masses. An additional resonance should therefore become visible in such energy 
independent analyses and the results in (83), (84) don't depend on the parameterization 
chosen for the fit. 

The present evidence for /o(1370) comes almost exclusively from fitting energy 
dependent parameterizations to the angular moments or Dalitz plots. For example, Bugg 
[125] has fitted the CM-I moments (and other data) to a superposition of resonances 
with a sizable /o(1370) contribution in tttt scattering. The tttt S wave amplitude shows 
an extra circle between 1000 and 1500 MeV. From the size of this circle in comparison 
to the elastic circle below 1000 MeV we conclude that the elastic width of /o(1370) in 
this fit is about Xtttt ~ 0.35 (figure 16 in [125]). Such an effect can hardly be consistent 
with the 7] — 6 "data" in figure 11 which in turn fit the original moments in all mass 
bins separately. Accordingly, although the overall fit to the moments (figures 14,15 in 
[125]) looks reasonable, there are significant local deviations between fit and data in all 
moments, also in the region around 1300 MeV. 

In the K matrix analysis by Anisovich and Sarantsev [37] /o(1370) appears with 
partial width x^.^^ ~ 0.19 — 0.26, depending on the solution (Table 6 in [37]), and it is 
larger or equal to the tttt width of /o(1500) with Xj^t^ ~ 0.19 — 0.23. Decay widths into 
TTTT of both states with equal strength are hard to reconcile with the energy-independent 
result (84) from CM-11 and the amplitude structure in figure 11. 

The GAMS Collaboration has performed a partial wave decomposition of their tt^tt" 
data [150] which has been used to reconstruct the S wave amplitudes (our figures 10 
and 11 from [102]). The Collaboration also determined the contributions of individual 
resonances to the S wave mass spectrum \S\'^ (figure 9). A model amplitude is fitted 
to the spectrum which represents a superposition of a real Gaussian- type background 
and Breit-Wigner resonances with arbitrary production phases. With the following 
production cross sections a x BR a good description has been obtained (see curves in 
figure 9) 

/o(980): 5.4±1.2nb, /o(1300): 70±15 nb, /o(1500): 12±3 nb, /o(2010): 3±1 nb (85) 

The fit requires the by far largest cross section for /o(1300) which is not at all visible 
in the amplitude analysis of our figure 10. This confiict arises using a fitting procedure 
which tries to derive a two-dimensional (complex) amplitude with several interfering 
components from a one- dimensional mass spectrum which ends up with ambiguities. 
One can anticipate that the model amplitude for the S wave would not fit the relative 
phases to the D wave represented by the /2(1270) resonant amplitude which enter the 
amplitude analysis. It should also be noted that the second dip near 1500 MeV in the S 
wave mass spectrum is generated not only by /o(1500) but also by the large background 
from the isotensor 5*2 (see figure 10). 
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1.2.6. Inelastic scattering nn — )■ KK, rjrj. Data on tttt — >■ KK have been analyzed as 
well. The experiments with highest statistics are : 

(i) Argonne experiment on TT~p — )■ K~K^n and vr+n — > K~K^p at 6 GeV for 
'f^KK < 1-55 GeV by Cohen et al. [92]. Their complete analysis yields 8 different 
solutions; they are all eliminated but one. They find one S wave resonance with 
mo = 1425 ± 15 MeV and width T = 160 ± 30 MeV. Prom a determination of the 
zeros of Det {K~^) they exclude the existence of more than one resonance below 1500 
MeV. A peak in the mass spectrum near 1300 MeV associated with a smooth phase 
behaviour is explained by the interference between the background from /o(980) and the 
resonance at 1425 MeV (figure 33 of [92]). Their Argand diagram suggests a negative 
relative coupling between background and resonance, as also concluded in [61] for the 
Brookhaven experiment. 

(ii) Brookhaven experiment on vr^p — ?> K^K^n at 23 GeV and for m^f^ < 2.4 GeV 
by Etkin et al. [91]. Their best fit includes a resonance at mo = 1463 ± 90 MeV and 
width r = 118±l}f in perfect agreement with /o(1500) from PDG. Their global fit 
with /o(980) and /o(1500) (but without /o(500)) represents the S wave phase 5q well, 
but not the peak near 1300 MeV. 

These results are summarized by Biittiker et al. [154] in figure 12. The phase shifts 
at the KK threshold join smoothly with those from elastic vrvr scattering which coincide 
below threshold because of unitarity (see, for example [155]). Near the KK threshold 
the presence of a P wave - the tail of the p meson - in tt+tt^ K^K~ (Cohen et al.) 
allows the accurate determination of the S wave. As there are conflicting results in 
this region a preference is given for low energies E < l.'i GeV to the data by Cohen et 
al. The tttt phase shifts obtained in a dispersion relation approach [100] also prefer the 
upper branch of phases which approach 225° at KK threshold. 

Subsequent descriptions in terms of resonances including the broad object /o(500) 
yield a good qualitative description of tttt and KK data without an additional resonant 
state near 1300 MeV: a K matrix fit [118] with three poles and a fit with the ansatz 
(83) [151]. On the other hand, a contribution from /o(1370) is suggested in [90] with 
an effect hardly visible in the phase motion. Indeed, the phase behaves very smoothly 
around 1370 MeV in the data (Left panel of figure 12), while there is a strong effect 
from /o(1500) similarly to the elastic tttt channel, and, of course, from /o(980). 

(iii) Finally, the inelastic process tt'^'tt" — )> rjrj has been measured by Binon et al. 
(GAMS Collaboration) [156] (more recent results in [157]). Also in this case, a significant 
effect is seen for /o(1500) with mass and width consistent with PDG. The state interferes 
with the broad background, opposite in sign to KK [61]. It forms a circle in the Dalitz 
plot also seen in [90]. A clearly visible signal from /o(1370) can hardly be identified. 

7.3. Decays of D and B into /q mesons 

In the weak decays of heavy quark mesons some well defined qq states evolve from the 
intermediate weak and strong interaction processes and they finally can form isoscalar 
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Figure 12. Phase and magnitude of / = 5* wave amplitude for tttt KK vs. E^r 
(data [101, 92, 91]); curves for Ett-k > 1 GeV are polynomial fits (figure from [154]). 
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Figure 13. tttt S wave amplitude and phase extracted from decays Df 7r+7r~7r+ 
(BaBar Collaboration [122]); right panel: Argand diagram for tttt amplitude (the phase 
is normalized to 7r/2 at the /o(980) peak). The points are not equidistant in mass. 



mesons. Then, these processes may help to clarify their intrinsic structure. 

7.3.1. Decay of Ds mesons. The dominant subprocess for /o production in the decay 
— )■ 7r+7r^7r+ is identified as 

Dj^n^ + ss; ss /o tt+tt". (86) 

Earlier studies of this process by E791 [158] and FOCUS [159] Collaborations have 
shown a strong signal from /o(980) and some resonant structure associated to /o(1370) 
or /o(1500). In a dedicated analysis by Klempt, Matveev and Sarantsev [160] of E791 
data no positive evidence for /o(1370) besides /o(1500) was found, although such a 
contribution has not been excluded. 

In a more recent high statistics investigation of the decay Dg — >■ Sti by the BaBar 
Collaboration [122] the known resonances p, /2 have been represented by Breit-Wigner 
forms while the 5* wave amplitude As is parameterized in a model independent way by 
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an interpolation between 30 points in the complex plane § 

Asim^n) = Interp(cfc(m^^)e^'^'=("*--))|fc=i...3o. 

The resulting amphtudes and phases are shown in figure 13. One observes clearly the 
peak from /o(980) and another peak below 1500 MeV. The phases show two regions of 
strong variation with mass, near 1000 MeV and near 1500 MeV (middle panel). Near 
1370 MeV the phase motion per mass interval is near minimum, without resonant signal. 

The Argand diagram (right panel) shows a large circle peaking at 980 MeV with a 
small circle near 1500 MeV superimposed curving outward. This behaviour corresponds 
to a negative relative coupling of /o(1500) with respect to the background from /o(980) 
(and some /o(500)) according to formula (46) with T ~ - Trc^^^b _ 

This is to be contrasted to elastic tttt scattering in figure 11 where the small circle at 
1500 MeV is curving inward corresponding to positive coupling. The same situation with 
different sign of coupling has been discussed before in tttt KK and tttt — )> rjrj [61] and 
can be interpreted in terms of the near flavour octet composition of /o(1500) as in route 
2 according to (34). There is no evidence for an extra circle in the Argand plot near 1370 
MeV. However, as the initial coupling to ss is not known one cannot derive a limit on 
the branching fraction of /o(1370) tttt. The shift of the /o(1500) peak towards smaller 
masses is due to the interference with the background amplitude from /o(980)//o(500). 
The negative relative couphng of /o(980) and /o(1500) to ss is also refiected by the clear 
depletion in the Dalitz plot of Dg — > Sn at the /o(980) — /o(1500) band crossing [122], 
and this has been seen already in the data from the E791 Collaboration [161], see also 
the discussion of phases in [70] . 

7.3.2. Decay — )■ K'^Jq, /q — )■ tttt, KK. The tttt and KK mass spectra from these 
decays as measured by Belle [82] are shown in figure 14. There is a clear signal of 
/o(1500) in the K^K~ spectrum, but not in 7r"'"7r~, despite the 4 times larger decay 
probability in favour of the tttt channel. Therefore, a new state /x(1500) has been 
suggested [82] which looks like /o(1500) except for the branching ratios. || 

This effect can also be explained [48] by taking into account the different quark 
composition of /o(1500) near fiavour octet as in (35) and /o(500) (or /o(980)) taken as 
flavour singlet (or glueball) along route 2 which leads to an opposite sign of interference 
in the strange and non-strange components: one can obtain a constructive interference 
of both components in K^K^ and then one finds a destructive one in tt^tt^. In this 
model represented by the curves in figure 14 the broad background from /o(500) and 

§ This method has been used first by the Fermilab E791 Collaboration [158]. 

II In the analysis of iJ — )■ KKK by BaBaR [162] the mass peak near 1500 MeV has been fitted by 
X(1550) but with a large width of 257 ± 33 MeV, twice as large as of /o(1500). We note however, that 
the fit - contrary to the analysis [48] - corresponds to an approximately destructive interference of this 
resonance with background. The data are not so well reproduced in the resonance region and another 
solution with smaller width may exist. 
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Figure 14. tt+tt^ and K^K^ mass spectra in B-decays (Belle [82]) in comparison 
with a model amplitude jTp of the coherent superposition of /o(980), /o(1500) and a 
glueball {gh) forming the broad background. Also shown are the individual resonance 
terms The background [gb) and /o(1500) interfere destructively in tttt and 

constructively in K~^K~ final states (figure from [48]). 



the narrow resonances are superimposed as in (46) following Michael [99] 

Tkk ~ + csT/^dsooe'''^" + C4Ty„(l5oo)e'*'«^ (87) 

Here Tj denotes Breit-Wigner amplitudes and 6gb the S wave phase shift corresponding 
to /o(500) which varies between 90° and 130° in the considered mass range. The 
phase difference of the KK component, i.e. of C2 and C4 in (87) near vr, explains the 
different structures around 1500 MeV: the strong peak at 1500 MeV in K^K~ and the 
small enhancement at 1400 MeV followed by a drop towards the resonance position at 
1500 MeV in the vrvr spectrum. The interpretation in terms of subprocesses has some 
complexity by the contribution of different diagrams [48] and we do not discuss this 
further here. 




Figure 15. Leading diagrams for decays ^ ^ J/^ V^v'^ /o,/o and i?" 



7.3.3. Decay 5° J/ipfo, /q — )■ tttt. This process has been measured recently by 
Belle [163] at the Y{5S) resonance in e~^e~ annihilation and with higher statistics by 
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Figure 16. Mass squared spectrum m^(7r+7r ) in i?° — J/V^Tr+Tr decays; the dotted 
and dashed hues indicate the background and signal contributions respectively, figure 
from [83] (LHCb CoUaboration) . 



the LHCb Collaboration in pp collisions at 7 TeV [83] . It proceeds dominantly through 
the intermediate steps 

Bs J/lp + SS] SS ^ fo ^ TTTT. (88) 

The dominance of the ss intermediate state has been explored in particular in connection 
with the rj — rj' mixing [164, 114, 165, 166] (see figure 15, middle panel). 

The TT+TT^ mass spectrum has been described by both Collaborations in a model 
fit with contributions from /o(980) and /o(1370) above a broad background in the nir 
S wave. The LHCb Collaboration also included /2(1270) in the minimal fit (at a 10% 
level) and investigated the inclusion of small contributions from other states which led 
to only minor improvements. The mass spectrum shown in figure 16 shows the narrow 
peak from /o(980) and an enhancement followed by a drop near the position of the 
/o(1500) resonance (m^^ = 2.25GeV^). In their preferred solution they find mass and 
width of the second resonance called "/o(1370)" as ^/^(isyo) = 1475.1 ± 6.3 MeV, and 
r/o(i370) = 113 ± 11 MeV. Comparing these numbers with the PDG listings we find 
the identification of this state with /o(1500) to be more appropriate: the width agrees 
well with rjQ(i5oo) = 109 ± 7 MeV while it is much smaller than the PDG estimate 
r/o(i37o) = 300 — 500 MeV. The mass is smaller by 2% in comparison with the PDG 
value m/Q(i5oo) = 1505 ± 6 MeV, but the mass may be affected by additional systematic 
effects from the interference with other contributions. Also, the "/o(1370)" mass is found 
at the upper edge of the large range suggested by PDG: mjQ(i37o) = 1200 — 1500 MeV. 
Therefore we identify the heavier resonance in tttt as /o(1500) and not as /o(1370) in the 
Bg decay. An additional experimental check is possible by comparison with the K^K~ 
channel using the respective known branching ratios; the proper inclusion of interference 
effects is important according to the lesson from B — > KKK above. Indeed, a signal 
at 1500 MeV is observed [167] but there is an overlap with /2(1525). We also note the 
similarity of the mass spectra in figures 13 and 16 with peak below 1500 MeV. 
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1.3.J^. Mixing o//o(1500) — /o(980) jrom Bg and decays. The branching ratios for 
Bg — 7> J/ipfo also allow tests concerning the intrinsic structure of these states. In route 
2 (sec table 6) they belong to the same multiplet and the coupling of ss into /o(980) and 
/o(1500) according to (35) is approximately like 2j\f^ and — 1/-\/3. This expectation is 
compared with the measurement of branching ratios in table 11. 



Table 11. Fit fractions of J/tf;fo{980) and J/V'/o(1500) for two solutions found 

by LHCb [83], also rescaled for unseen decays (using B(/o(1500) nn) = 0.349±0.023 
[4] and B(/o(980) -> tttt) = 0.75±0.12 [169]; predictions (route 2) without phase space 



corrections; (/o//o ^ 


((- sin (j)sc)/ cos (l)sc 










preferred solution 


7r+7r~ fraction % 


corrected % 


route 2 


Phase 


B° ^ J/^/o(980) 


107.1 ±3.5 


214.2 ±34.9 


c.(. 


^273)' 


0. 


^ J/V'/o(1500) 


32.6 ±4.1 


141.6 ±20.1 


:^(- 




(241.5 ±6.3)° 


NonResonant 


12.84 ±2.32 


19.26 ±3.48 






(217.0 ±3.7)° 






0.(i()± 0.14 


O.O 


X PS 




Alternate solution 


5° ^ J/V/o(980) 


100.8 ±2.9 


201.6 ±32.5 


^(. 


v/273)2 


0. 


^ J/V/o(1500) 


7.0 ±0.9 


30.4 ±4.4 


^(- 




(181.7 ±8.4)° 


NonResonant 


13.8 ±2.3 


20.70 ± 3.48 






(232.2 ±3.7)° 


/o(1500)//o(980) 




0.15 ±0.03 


0.5 


X PS 





The Belle Collaboration finds the branching ratios B{Bg — ?> J /ipf o{980)) = 
l.lQloil X 10"^ and 5(5° ^ J/^/o(1370)) = 0.34+°;}^ x 10"^ if we add the errors 
in quadrature. After corrections for unseen decay modes we find (taking the second 
resonance as /q = /o(1500)) 



p35(fiO J/v:'/o(1500)) , ■ 
^/i//o = (p/)3ij(^o^ j/^/p(980)) ' ""/o/Zo = S 0-15 ± 0.05) LHCb alternate solution, 



1.18 ± 0.25) LHCb preferred solution, 
0.15 ± 0.05) LHCb 
1.12 ±0.61 Belle. 



(89) 

with p{p') the momenta of /o(980) (/o(1500)) in the Bg system. Combining the LHCb 
preferred solution with Belle yields 



f^/fo 



1.17 ±0.23. (90) 



It is interesting to compare this result with the corresponding ratio for decays into 1] and 
f]' which should agree approximately for route 2. This ratio has been measured recently 
by [168] as i?^ = B{Bg J/^r]')/B{Bg -> J/^prj) = 0.73 ± 0.14 (stat.) ± 0.02 (syst.) 
which leads to the phase space corrected ratio 



p;, B{Bg ^ j/xbf]) 

pI B{Bg ^ J/ijT]') 



K/,' = 1 o. o^ . = 1-10 ± 0-22. (91) 
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This number is in a remarkable agreement indeed with the result in (90) supporting the 
similar flavour mixing of /o(1500) — /o(980) and 77 — 77' with inverse mass ordering as in 
(35) for route 2. 

These ratios can be compared also directly with the expectations according to the 
quark-gluon structure of these mesons. If we allow for a gluonic component with mixing 
as in (18) we find (see also [166]) 



tan^ - 



'ps 



COS 



2 



COS (/)^ = 0.86 ±0.09. (92) 



where (f)ps = 42° has been used. This is easily compatible with having no gluonic 
component in i]', alternatively, one can admit a contribution with (pQ ~ (30.7^^2^5)°. 
Accordingly, introducing a gluonic component for /q = /o(1500) as 

1/0) = cos^gKcos^scI?^^) - sin0sc|ss)) + sincpalgg) (93) 
we find for route 2 

r^,/^^ = cos^ tan^ 0^, ^ cos 0§ = 1.87 ± 0.32. (94) 

where we have used our result 0sc = 30° in (68) and r^'/f^^ = 1-17 in (90). We find that 

f " 
this result for cos (J)q exceeds the physical region by 2.7a. This discrepancy may indicate 

a deficiency of the model with gluons mixed only into /o(1500) (if not of the /q — /q 

mixing model altogether), but looking at the success of the relation r^^^.^^, ^ ''^^'/fo '^^^ 

could also think of an underestimation of errors in the branching fractions of Bg decays, 

of the correction factors for unseen decays or our mixing angle (psc- An interesting test of 

the resonance analysis with /o(1500) would be the corresponding analysis of the decays 

Bg — > J/ipKK where different interference phenomena appear, as we could witness 

already in the B ^ K + X decays (see figure 14). 

The production phase of /o(1500) measured with respect to /o(980) in the ss 

component selected here is expected in route 2 as = 180° which is near the observed 

results in table 11. For illustration, consider our approximate representation (35): in 

the nn — ss coordinate system /o(980) = (1, \/2)/\/3 and /o(1500) = (\/2, —l)/\/3, so 

the relative phase for the ss component is = tt and for the nn (or dd) component 

it is A^'' = 0. The phase in the LHCb preferred solution is found as 

A(j)' = (241.0 ± 6.3)° (95) 

with an additional shift by ~ 60° which could be related to a moving background phase 
as in the model of the type (87) (the /o(500) component in figure 11 moves from about 
80° to 140°). 

The flavour mixing /o(980) — /o(1500) can be measured also from Bd decays, a 
measurement which may become feasible in the near future. In close analogy to the 
Cabibbo-suppressed decays B^ — >■ J/iptj^''^ [164] one can study the corresponding decays 
into /o//o- The dominant subprocesses in this case are (figure 15, left panel) 

Bd J/^p + dd, dd fP, r]^'\ (96) 
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One obtains for decays in the standard mixing scheme 

p3, BjB, J/i^rj) _ 
'^1^' - v% B{B, ^ JIH) ' " ^ ' 

In the same way, assuming an /o(980) — /o(1500) mixing as in route 2 

, _ p^BjB^, ^ J/^/o(1500)) , _ , 
"^fo/fo - J/V./o(980))' "^^o/Zo - '^^ ^^^) 

with the same mixing angle (psc- For these mixing schemes one obtains the constraints 

Also in route 2 one expects (pps ~ or 

4//o^</.'- (100) 

For the production phase of /o(1500) relative to /o(980) we obtain in B^ — )■ J/i/^nn 
decays A^*^ = 0, or, above a background with moving phase as in tttt scattering 

A(f)'^ ^ 60°. (101) 

7.5.5. Mixing angles of isoscalars /o from B^ and Bg decays. The determination of the 
individual flavour mixing angles of isoscalars like /o(980) and /o(1500), but also of 
/o(500), i.e. the ratio of nn and ss amplitudes of these states, is possible from the ratio 
of Bg and decays. As in case of ri,ri' [164], one can obtain the mixing angle (pgc for 
/o(980) or the mixing angle 0'^,^ for /o(1500) from the ratios 
p'gTjB, ^ J/^/o(980)) 1 

plTjB, ^ J/^/o(1500)) 1 

= p3r(E3 ^ j/^/o(i50o)) = 2^^^--) • 

where Fqkm — Vcd/Vcs is calculated for the dominant tree diagram and penguin diagram 
contributions and with the CKM matrix elements from PDG one finds \Fckm\ — 0.225. 
The mixing angle (f)sc for /o(980) should agree with our determination in (68), the angle 
(pg^ for /o(1500) may be different if the two states arc in different multiplets. 

Under the assumption that /o(1500) and A'*'^(1430) belong to the same multiplet, 
the mixing angle for /o(1500) can be determined also from the ratio 

_ Pi ^ J/^/o(1500)) _ 

" Pinsd ^ j/vi^*°(i43o)) - (^°^^ 

in analogy to a corresponding ratio of r{Bs — > J /'4>rf) /T{Bd — >■ J/'4^K^) (see figure 15). 
It is interesting to note, that these determinations of mixing angle between the different 
flavour components are independent of a possible presence of a gluonic component. 

In the same way one could also study the flavour content of /o(500)/(T. 

We have no direct predictions from the other classification routes discussed. Some 
aspects of the tetra-quark model for /o(980) are elaborated in [114]. 
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14. DecayofJ/i^^<t>{lQ2Q>)fo, /o ^ tttt, KK 

In the TTTT final state of this J/-^ decay the BES Collaboration [116] observed a peak 
around 1335 MeV followed by a drop near 1500 MeV. In their model fits they include, 
besides /2(1270), the scalar /o(1500) and require /o(1370) as well with mass m — 1350±5 
MeV and width T = 265 ± 40 MeV. 

The J/ip decay is expected to proceed dominantly (according to the "singly 
connected diagram") through a gluonic intermediate state and, triggered by 0(1020), 
the vrvr final state evolves from an intermediate ss configuration. Therefore, we expect 
a similar final state interaction as in the decay Bg J/tlJini (see discussion in section 
7.3.3 and figure (16)) and in tttt — > KK (see discussion in section 7.2.6 and figure 12). 
Indeed, the tttt spectra in these processes with an enhancement near 1400 MeV and a 
drop near 1500 MeV look all rather similar. There, evidence was only given for a single 
resonance interfering with some smooth background amplitude and the resonance was 
compatible with /o(1500). It is therefore important for a definite conclusion about an 
extra state "/o(1370)" to analyze the observed peak in J/i/j ^ (p+mr decays with higher 
statistics in a bin-by-bin phase shift analysis. 

In the K^K~ channel the BES data [116] do not show any enhancement around 
1350 MeV and the branching fraction ratio 

- S(/o(1370) ^ tttt) - ^ 

has been determined. Other measurements of this ratio hsted by PDG [4] are spread 
over a range up to Rkk/^-k i$ 1-0) therefore 

5(/o(1370) ^ KK) < 0.1. (106) 

The measurement r(7777)/r(47r) = 4.7 ±2.0 x 10~^ [170, 4] implies this number to be an 
upper hmit for S(/o(1370) — > 7777). So, as an estimate 

5(/o(1370) ^ 2 pseudoscalars) < 20% (107) 

and almost all decays should go into higher multiplicities. 



7. 5. Central production in hadron hadron collisions 

Central production processes are selected by requiring two gaps in rapidity y — ^In ^3^, 
empty of particles, adjacent to the direction of the incoming hadrons, and smah 
momentum transfers t — {pi—pfY between the incoming and outgoing forward hadrons 
a and h (see figure 17). For such processes a description in terms of Regge theory 
applies. The exchange of Reggeons IR or Pomerons IP determines the dependence of 
the respective cross sections on the primary energy ^/s. At high energy ^/s the Double 
Pomeron Exchange (IPIP) process dominates with an approximately constant cross 
section whereas subleading contributions fall like a{IPIR) ~ Vv^ and a{IRIR) ~ 1/s. 
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a a 
[R,rP 

[R.rp 

b b 

Figure 17. Central production of meson M in a + 6 collisions with two rapidity gaps 
and exchange of Reggeons and/or Pomerons. 

While the Reggeon is considered a qq object, the Pomeron is thought rather to be 
dominantly gluonic, more discussion follows below in section 8.1. Therefore the strength 
of central meson production can give a hint towards its possible gluonic component 
[45, 46]. In particular, a search for scalar mesons in many different channels has been 
performed. 

7.5.1. The nn system, has been studied extensively at the ISR (Axial Field 
Spectrometer) [171], at the SPS coUider by WA102 [172] and by GAMS [173], and 
now at the LHC by ALICE [174]. In figure 18 we show as an example an amplitude 
analysis of the tt^tt'^ system including S and D waves which has been carried out using 
experimental moments (Y/^) and relations (80). Two distinct solutions appear. Solution 
2 with a D wave at threshold much stronger than S wave is considered unphysical. 
The following features of the physical solution 1 in figure (18) are noteworthy: 

(i) Low mass peaks in S and D waves. In the S wave intensity there is a strong peak at 
low mass m^Tr ~ 400 MeV but without associated phase variation. In addition, there is 
also a peak in the Dq wave at 500 MeV much below the first known resonance /2(1270), 
very different to elastic tttt scattering, for example. As pointed out in [175] a very similar 
phenomenon has actually been observed in 77 7t~^7i~ (not so in tt^tt^), where the low 
mass S and D wave peaks are related to the dominant one pion exchange process at 
low energy with direct coupling of photons to charged pions. The pion exchange pole 
near the physical region of scattering at t = causes the enhanced forward-backward 
scattering generating the large D wave at first. This analogy suggests the appearance 
of one pion exchange in the subprocess 

IPIP TTTT 

as well, which adds to the direct resonance formation as in figure 17. A Regge calculation 
of this exclusive process has been presented recently [176]. The presence of the broad 
/o(500)/(T is signaled by the interference with /o(980) (sharp drop in |5'op at m^T^ ~ 1 
GeV) as in tttt elastic scattering. At higher energies other exchange processes like p- 
exchange are possible [5] . 

(ii) There is no phase coherence. Contrary to the reaction n~p — )■ 7r°7r°n [150] the phase 
differences ^pso ~ ^^"^ V^^o ~ Vd^ very different. In particular, while \Dq\'^ and 
\D^\'^ both show a clear signal from /2(1270), only tpsa — ^r,- I'eAects a related phase 
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Figure 18. Partial wave amplitudes for the centrally produced n w system in 
pp p{n°n°)p at = 450 GeV, figure from [172](WA102 Collaboration). 



variation from its interference with the slowly moving S wave, whereas (pso — is 
without resonance structure and so "phase coherence" is strongly violated. 

Therefore, one cannot construct a model based on formulae (80) with a resonance 
/2(1270) in Dq and interfering with the same -S"^. The violation of phase coherence 
as in figure 18 is expected if /2(1270) couples to different partial waves L^^^ with different 
nuclcon helicities. As an example, let the Sq wave couple only to Sq^_^_ ^_|_ but let 
/2(1270) couple to Dq"^^ and also to D^^_^_ (the lower (+) and (— ) refer to the 
nuclcon helicities), then /2(1270) appears in both l-^fP moduli, Sq-Dq is resonant 

while So ■ Di = similar to observations. 

Such effects can occur if several exchange mechanisms are at work with different 
spin composition and "spin coherence" is violated. Then, in (80) it may be necessary 
to replace the single terms 

\S\^ ^\S\^ + \S\^; 

Re{SD-*) Re{SD-*) + Re{SD-*) (108) 

where the terms S and S correspond to contributions with different nucleon helicities. 
Obviously, by addition of arbitrary incoherent contributions as in (108) one cannot in 
general solve for all amplitudes in an energy independent analysis, but possibly a model 
with selected contributions to describe all moments can be constructed. 
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In a more general approach one may include a description for the transverse 
momentum dependences of the helicity amplitudes and in that way provide additional 
information for a full separation of partial wave amplitudes. Transverse momentum 
correlations of outgoing protons have been suggested as tool to discriminate different 
intrinsic structures of the centrally produced object ( "glueball-filter" ) by Close and 
Kirk [177]. Different effects have been found indeed for the resonances in the tttt, KK 
data obtained by the WA102 Collaboration [178]. Especially, the azimuthal final pp 
correlations behave differently for /2(1270) and /o(1500) [179] which implies again 
that a model with "spin coherence" as incorporated in the amplitude analysis above 
is incomplete. The t-dependcncc of helicity amplitudes in the general framework of 
Regge theory and for specific models has been discussed by Kaidalov et al. [180]; and 
so this framework may be used for an extended analysis including helicity dependent 
^-distributions. 

(iii) Contributions from states /o(1370) and /o(1500) have to be determined from 
a partial wave analysis, i.e. from the complex amplitudes as discussed in subsection 
7.2.5. In view of the previous remarks this requires an extended analysis not relying 
on spin-coherence. Possibly at the higher LHC energies spin coherence is again a good 
approximation. In any case, a determination of the resonance fractions from |5'o]^ alone 
without inclusion of phase movements does not yield reliable results as was demonstrated 
in subsection 7.2.5. For this reason the results from WA102 on /o(1370) 7nT,KK 
[172, 178] are not considered here as evidence for the existence of this state. 

7.5.2. Production of /o(1370) and /o(1500) in the in system. The An decays of 
the scalar mesons have been investigated extensively in central diffractive production 
pp -)■ p{An)p (WA102 [181, 182]) and in pp 57r° and pn n~An° annihilations 
(Crystal Barrel [183]). A comparative discussion is given by Klempt and Zaitsev [5]. 
They have collected the decay branching ratios of both resonances into different An 
components as measured by central production and by pp annihilation experiments (see 
table 12) 

Table 12. Decay branching ratios of isoscalar mesons into An channels from central 
production (CP) and pp/pn annihilation (from [5]). 
# taken from PDG listing 





/o(1370) 

pp/pii CP 


/o(1500) 

pp/pii 


CP 


r47r/rtot 


0.80 ± 0.05 


0.73±0.21# 


0.48 ± 0.05 


Tpp /r47r 


0.26 ±0.07 ~0.9 


0.13 ±0.08 


0.74 ±0.03 




0.51 ±0.09 ~0 


0.26 ±0.07 


0.26 ±0.03 



Some discrepancy for /o(1500) disappears, if the different decay channels included 
in the analysis are taken into account. In case of /o(1370) essential discrepancies remain. 
Especially, in central production there is a clear signal from /o(1500) in the aa channels 
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(including 47r°) but no obvious signal from /o(1370) (cr stands here for low mass tttt S 
wave). On the other hand, in pp {pn) there is a strong peak in the mass spectrum 
at ~ 1470 MeV but with a larger width than expected for /o(1500) and this effect has 
been related to /o(1370). The analysis has to take into account the contributions from 
different processes and decay chains in Stt processes and combinatorial backgrounds in 
57r° not present in central production. An attempt to analyze the phase motion of the aa 
decay has lead to the conclusion [5] that there is only one resonant state with a motion 
consistent with /o(1500); the strong pp enhancement at threshold in central production 
which has been related to /o(1370) and its complete absence in aa in this view could 
be a consequence of p-Regge exchange in IPIP scattering, similar to the threshold 
enhancement in tttt which has been attributed to pion exchange as discussed before. 
This mechanism would explain the absence of /o(1370) — )■ aa in central production and 
similarly the absence of the rji] production where a suitable Regge exchange is lacking. 
There is also the interesting phenomenon of a dip in the pp mass spectrum which moves 
depending on the kinematic selection for the outgoing protons [181]. For a particular 
selection (azimuthal angle between protons 135° < < 180°) the deep sharp dip is just 
at the mass of 1500 MeV which suggests a negative interference of /o(1500) with a broad 
background reminiscent of the similar phenomenon in tttt scattering (figure 9) . 

While the signature of /o(1500) is very clear in (47r), especially in (aa), there is no 
compelling evidence of additional contributions which would require a new resonance 
at lower mass around 1370 MeV. As learned in the 2 particle channels, instead of the 
broad /o(1370) (F > 300 MeV) there may also be a contribution from the long tail of 
/o(500 - 1000) with F ~ 1000 MeV. 

7.6. Flavour structure o//o(1500) from its decays 

This state has been considered at first as a glueball candidate following a prediction 
from quenched lattice QCD. The study of decay branching ratios suggested a mixed 
state of qq (mainly octet) and gg as discussed in route 1 above [35, 51]. A pure qq 
composition was found incompatible with certain decay properties [184]. On the other 
hand, the hint towards a negative phase between K^K" and t:^t:~ decay amplitudes 
[61, 118] limits the contribution from a glueball component which would prefer a singlet 
type qq state. 

A direct measurement of the mixing angle a (= (psc in route 2) independent of the 
gluonic component is possible through the measurement of decays (see subsection 
7.3.5). The negative relative sign in D,, decay between ss — )■ /o(980) — )■ tttt and 
ss /o(1500) — )■ TTTT (see subsection 7.3.1) requires an octet type mixing angle a = (psc 
in (32) with 

tan0,c>O (109) 

once the singlet type mixing for /o(980) is established. In the following we will 
investigate in more detail the consequences of the measured decay branching ratios 
for the constituent structure of /o(1500). 
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The decay couplings jfj of /o(1500) depend on the qq and gg components through 
the mixing angle a (see (15), (32)) and the parameter rc ~ G/g, as shown in table 3. 
Now we consider ratios Ri of branching ratios Bk as in [35] . From the PDG fits [4] using 
the phase space correction 7^ = B^/qk we derive the ratios 

R^^jWm^VVi ^ 0.208 ±0.039 (110) 

7'[/o(1500) ^tttt] ^ ^ 

72[/o(1500) ^tttt] 

72[/o(1500) ^ , , 

^3 = o\ r\..r.L r = 0.321 ± 0.034 112 

72[/o(1500) ^ tttt] ^ ^ 

(using also g^^/ = 195 MeV for R2 [35]). For a given experimental ratio Ri we can insert 
in (110)-(112) the results for 7? from table 3. Then, for each value of the mixing angle 
a there are two corresponding solutions for the glue contribution Vq-. Solving for Vq- 
we find (with = and p — R — 1) 

Vq-^ {-y/SRi^ 1) = —\^3Ri cos a ± (cos^^cosa — -\/2sin^ sina) (113) 
V 3-^2 = — •\/3i?2 cos a ± -\/2 sin cos 0(cos a + -\/2 sin a) (114) 
r^3 (\/3i?3 ^ 2) = — A/3i?3 coscK ± (cos a — v^sinct) (115) 

In figure 19 (left panel) these solutions r^- for every R^ are shown for the given a. There 
are two special points for tana = —l/y/2 where all solutions coincide in a trivial way if 
the numerator Ni and the denominator for RiDi — Ni vanish: 

ai ^ -35°, a2 ^ 145°; with rai ^ -0.82, rG2 ~ 0.82. (116) 

Otherwise, there is no crossing of the solutions from the three ratios Ri expected for an 
exact solution for ra, but rather some pairs of ratios come close to each other {rQ2, Tq^) 
and (r^^r-g). 

In the right panel of figure 19 we show, how the varies with the glue component 
Tg for a mixing angle 30°; these two parameters determine the ratios Ri according to 
table 3. Also shown are the contributions to from Ri separately which vanish at the 
exact solutions for r^^ seen in figure 19 (left). There is a minimum 

forces 30°: r^^ +0.7 ±0.1. (117) 

For this solution there is a good agreement with R2, R3 but a discrepancy with the rjr) 
channel (i?i) by > 2a. There is another solution for ro ~ —0.4 but with an unacceptable 
> 30. The minimum is found to be practically independent of a if ra moves along 
the solutions (r^a^'^Gs)- that path there is also the solution (a = 0°, rc ~ 0) 
for /o(1500) without glue. Therefore, it needs further information to determine the 
mixing angle a uniquely. Some information comes from the 27 decay rate (see, for 
example [184]) which needs some model assumptions (see next subsection). If /o(980) 
and /o(1500) belong to the same multiplet as in route 2 then there are further relations 
as well as discussed above. A direct determination of the mixing angle is possible by 
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Figure 19. Glue contribution ra in the /o(1500) state as function of mixing angle 
a — (f)sc for fixed ratios Ri of branching ratios as in (110)-(112) (left); vs. ra for 
observed ratios Ri and full at mixing angle a = 30° (right). 



comparing i?^ and decays as in (103). Further information is available from the 
interference with other states like /o(980) or /o(500 — 1000) in different channels. 

According to Anisovich [185] the decay coupling go of the qq state into mesons 
and the coupling Gq for the gluonium decay occur in the same leading order of a l/Nc 
expansion. Then the couplings in (24) are of the same order 

90 ^ Go. (118) 

For the allowed solution of following rQ2, ^g3 deduce from the figure 19 the bounds 
~1 ~ '^G ^ 1 and with condition (109) from the /o(980) — /o(1500) interference 

0<rG<l (119) 

so for /o(1500) the gluonium contribution ~ G does not exceed the quarkonium 
contribution ^ g. As an example, if we assume the mixing a = (psc ~ 30° as in 
route 2 (^^^ ~ 4>ps) then we find the glue component in (117) 

G/g ^ tan 0g ~ 0.7; 0g ~ 35°. (120) 

The alternative solutions Tq^jTq^ have typically opposite sign and are disfavoured by 

7. 7. Production in 77 scattering 

The Belle Collaboration [94] has measured the angular distribution of the process 

77 — > 7r°7r° 

with high statistics. A partial wave analysis of the 7r°7r° system in the cms energy 
range 0.6 GeV < W^^ < 1.6 GeV has been carried out with S, Dq and D2 waves 
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included. As the photons arc Tinpolarized some ambiguities remain. A parameterization 
in terms of resonances and backgrounds yields parameters for /2(1270) and in the S 
wave for /o(980) and an additional resonance Y, either /o(1370) or /o(1500). They 
obtain the mass m = 14701^1^55 MeV and the total width T = 90lil^^ MeV and partial 
width T^^B^foiY) 7r°7r°) = lllglr^^ eV. The large systematic errors are related to 
uncertainties in the coupling of helicity pL — oi /2(1270) and the KK/nn branching 
ratio of /o(980). Mass and width are nicely compatible with PDG values for /o(1500) 
while for the alternative /o(1370) the width is smaller than the PDG estimate 300 — 500 
MeV. 

With the identification of the extra state as /o(1500) we can correct the above 
radiative width for other hadronic decay modes and obtain 

/o(1500) : r,,/r,,, = 95tim' eV, T!^f'''^/T;f''^ = 0.38l°:i^,l^:°« (121) 

We may also compare with constituent models. Assuming a mixing as in rj according 
to route 2 with /o(1500) ~ (uu + dd — ss)/y/3 we obtain the ratio, following [61] 

^°°)/r-(9«°) = i X {rufjma^r eV 

= 0.53 / 1.8 / 4.0 eV for p = -3 /I /3 (122) 

There is some uncertainty in the extrapolation in mass from the ao(980) to /o(1500). In 
Born approximation which is good for light pseudoscalars the power in (122) is p = 3; 
such a power is also consistent with the radiative widths of the lightest tensor meson 
nonet if the mixing angle ccr = 81° is used [3]. Taking only phase space with p ~ 1 
or a form factor with p = —3 wc find the lower values in (122). Given the large 
upper systematic experimental error there is no real contradiction with the pure quark 
composition models, on the other hand, there is enough room for gluonic contributions 
which would lower the qq predictions accordingly. 



8. Gluon rich processes in comparison 

So far, the different routes for classification of scalar mesons have not yet led to a 
generally accepted identification of a light glueball. In route 1 and 3 there is the need 
for /o(1370) where a clear evidence is missing: there is no established 2-body decay 
and their total decay rates are strongly limited (10-20%); the decay rates into higher 
multiplicities are controversial. A good argument why 2-body decays should be strongly 
suppressed in comparison to /o(1500), for example, is lacking as well. New experiments, 
for example GOMPASS [186] could reach higher sensitivities in such studies. 

In route 2, the asymmetry between a and k is a bit puzzling, although not fatal; 
/o(500) may not be purely gluonic and the mixing properties of isoscalars have to be 
better understood, possibly from Bd, Bg decays. 

The existence of glueballs as a basic prediction of QCD is now dependent on the 
existence of a supernumerous state in the classification which is found controversial. 
Therefore there was the hope to get an additional signal from the enhanced production 
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of glueballs in "gluon enriched" reactions. Looking over the discussions of the previous 
sections, there is not such a really striking predominance of a particular isoscalar state 
in the gluon rich processes (i)-(iii) listed in subsection (3.2.1). 

The reason could be that the relevant energies are of the order of the glueball mass 
in question itself in the range 1-2 GeV which may be too low for a gluon to appear as a 
well distinguished object. Rather, because of the running strong coupling the presence 
of intermediate qq components at low energies may be important as well. As an example, 
we consider the double Pomeron process (ii) and then we investigate processes where 
the presence of a gluon as primary agent is better established. 

8.1. The Pomeron as an example for mixed gluon and quark components 

At LHC energies of 7 TeV the double Pomeron process should dominate over Regge 
exchange processes for sufficiently large rapidity gaps. The Pomeron has isospin J = 0. 
If it is primarily a gluonic object we expect the emergence of gluonic mesons in such 
processes. A first result from the ALICE Collaboration with gap selections is shown in 
figure 20. While the vrvr spectrum without gap selection shows the p meson as dominant 



x1 0'' 




M(7nc) (GeV/c^) 

Figure 20. Central production of pion pairs at the LHC at 7 TeV with and without 
gaps in pseudorapidity —3.7 < r/ < —0.9 and 0.9 < i] < 5.1 (not corrected for 
acceptance), figure from [174] (ALICE Collaboration). 



structure, this state disappears if double gap events are selected. This is expected for 
the double Pomeron exchange process which cannot yield a final p meson with isospin 
1 = 1. Now, for the double gap events, the previously small structures related to /o(980) 
and /2(1270) become dominant. Certainly, the /2(1270) is a "classical" qq state as is 
nicely confirmed by the 2-photon decay rates according to the quark model (see, for 
example, [3]). The /o(980) is not typically considered as gluonic, though a small gluonic 
component cannot be excluded (see section (6)). Therefore, this result confirms that 
the double Pomeron process selects isoscalar states, but there is no evidence for the 
selection of gluonic states in particular. Such a behaviour may be related to a Pomeron, 
which is constructed from mixed glueball and qq trajectories as it is realized within a 
non-perturbative QCD approach [187]. 
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In this context, it is interesting to recall the results from the studies of deep inelastic 
scattering at the HERA ep coUider by the HI [188] and ZEUS [189] Collaborations on the 
partonic structure of the Pomeron as obtained from events with a rapidity gap adjacent 
to the target proton. At large momentum fraction z of partons inside the Pomeron 
the valence component is expected to dominate. It is found, that for z > 0.8 (for all 
scales presented with ^ 6 GeV^) the flavour singlet quark density dominates over 
the density of gluons. In a perturbative QCD treatment such dominant quark densities 
are initiated by Regge-like contributions and evolve according to the QCD equations 
for the quark and gluon densities [190] . Only in a global view, the gluonic component 
dominates: the integrated gluon density is much bigger than the quark density and the 
gluons carry a larger fraction of the total momentum as compared to the quarks. 

8.2. Processes with identified gluons 

In the gluon-rich processes (i)-(iv) the existence of gluons as intermediate partons is 
required by theoretical arguments but at the low energy we cannot be sure about the 
competing role of quark processes. We are therefore looking after processes involving 
higher gluon energies where the primary identity of gluons is better known. 

This is the case in the high energy production of quarks and gluons which evolve 
into hadronic jets. The identity of quark jets became obvious in e'^e~ annihilations 
above around ^/s ~ 7 GeV at SPEAR [191] whereas gluon jets have been identified 
around 30 GeV at PETRA [192]. With increasing jet energies, the quark and gluon 
jets can be better distinguished (see, for example [193], for a review). At LHC the 
knowledge of parton distribution functions allows the calculation of rates for quark and 
gluon jets and the selection of jets with a high purity (^ 80%) of either quark or gluon 
jets [194]. Then the fragmentation into particles of high momentum fraction z may yield 
information about gluonic mesons (sec below in section 9). 

Of special interest are exclusive 2-body channels (z ~ 1) which can be more easily 
identified also at lower energy and may carry a primary quark or gluon. An example is 
the weak ("penguin") decay of the heavy h quark (iv) b ^ sg which may contribute to 
B — > Kfo — > Knn. This processes have been studied in [48] . There are interesting tests 
for the members of the scalar meson multiplets but the isolation of the gluonic process 
has not been achieved so far. 

A promising study concerns the P wave charmonium (or bottomium) Xc (or xt) 
decays (v) which proceed through x ^ 99 with primary gluons (or other C = +1 heavy 
quarkonia) to which we turn in section 10. 

9. Leading systems in gluon jets 

A promising possibility to identify gluonic mesons is provided by the comparative study 
of leading particle systems in quark and gluon jets (see [195] for a recent discussion). 
According to the well known concept of quark fragmentation the leading particles at 
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large momentum fraction "Feynman x" are those which carry the primary quark of the 
jet as a valence quark, see figure 21. For example, leading particles in a m quark jet are 




Figure 21. Leading qq and gluonic particles in quark and gluon jets. 

n^{ud) or tt^{{uu + dd}/\/2) with half strength whereas 7r~ or K~ without valence u 
quark are suppressed at large x. In analogy, one can consider the fragmentation of a 
gluon and hypothesize that the leading particle has a gluonic valence component. 

Models of this kind with leading glueballs in jets, and also with leading isoscalars 
like ri,ri',u at large x have been suggested already long ago [196], x- distributions have 
been considered in [197]. The model for isoscalar qq meson production has been studied 
at LEP with quark and gluon jets [198], but no clear experimental support for the model 
has been established [195]. This means that the isoscalars rj and t]' behave like other qq 
mesons in fragmentation. 

A critical test for the presence of gluonic objects in the leading system of a jet 
is the measurement of charge and mass of the leading cluster Qieading and Mieading 
beyond a rapidity gap. Their distribution is expected to reflect the mechanism for 
colour neutralization [199]. In a non-perturbative approach an initial pair of qq with 
colour triplet charges creates another qq pair to neutralism the colour charges beyond 
the confinement radius Rc- For an initial gg pair the colour may be neutralized either by 
another gg pair ("colour octet neutralization") or by creation of two qq pairs ("colour 
triplet neutralization"), see figure 22. In a high energy jet one expects the charge 
distribution of the leading cluster to approach with increasing rapidity gap a limiting 
behaviour with charge Qieading = for colour octet {gg) and charge Qieading = 0, ±1 
for color triplet {qqg) neutralization according to the minimal configuration of partons 
beyond the gap. The mass distribution of the leading system with Qieading = should 
signal the presence of glueballs if any. 

9.1. Rate of neutral clusters at LEP above Monte Carlo expectations 

The distributions of charge and mass for the leading clusters in quark and gluon jets 
have been measured in e+e" annihilation at LEP by the DELPHI [200], OPAL [201] 
and ALEPH [202] Collaborations. The results by DELPHI on the distribution of the 
leading charge beyond a rapidity interval Ay = 1.5 are presented in figure 23. For the 
most energetic jet in a 3-jet event ("jetl", in 90% of quark jet) the distribution 

of the leading charge (sum of charges "SQ") follows closely the Monte Carlo expectation 
(JETSET [203]). This numerical calculation is based on the QCD approach to parton 
cascade evolution and string hadronization after tuning free parameters by experimental 
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Figure 22. Colour neutralization in quark and gluon jets (figure from [199]). 



data. On the other hand, for the jet of lowest energy ("jet3", a gluon jet with purity 
70%) there is a significant excess of events with Qieading = 0, clearly visible in the lower 
figures for the difference "data-MC". This excess is also shown in [200] to be built up 
by increasing the rapidity gap Ay. 
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Figure 23. Distribution of the sum of charges "SQ" (" leading charge") beyond a 
rapidity gap of Ay = 1.5 in gluon dominated (jet3) and quark dominated (jetl) jets in 
comparison with Monte Carlo (JETSET); analysis by DELPHI, showing the surplus 
of SQ = events in gluon jets (figure from [200] ) . 



In the analysis by ALEPH [202] gluon jets of high purity have been selected from 
the 3-jet events bbg with identified 6-quark jets. The distribution of charge has been 
studied in gluon jets without gap and in gluon jets with a rapidity gap. As can be 
observed in figure 24, the charge distribution in the gluon jet without gap follows the 
MC calculation (JETSET) whereas there is a large excess by about 40% in the jet with 
a rapidity gap. Models with non-standard colour connections along the parton cascade 
fail to reproduce the data as well. 

Finally, the distribution of masses of the leading cluster in figure 25, left side, shows 
again a good agreement with the predictions from the JETSET Monte Carlo for quark 
jets (right column), whereas for gluon jets the excess of events with leading charge 
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Figure 24. Left: distribution of the leading charge Qieading in a gluon jet of high 
purity without gap; right: the same for gluon jets with rapidity gap Ay = 1.5 showing 
the surplus of events for Qieading = over the JETSET Monte Carlo; other models 
(GAL, AR) with "colour reconnection" , figure from [202] (ALEPH Collaboration). 



Qieading = is observed at small masses, mainly for M < 2 GeV. An excess of events at 
small masses is also observed by OPAL (figure 25, right hand side). The mass spectra for 
selected charges with Qieading = don't indicate striking resonant structures as expected 

for narrow glueballs, a small enhancement around 1500 MeV in the (H 1 — ) spectrum 

is indicated while for (H — ) the mass distribution is smooth. Interestingly, there is no p 
signal which the models seems to predict; this could indicate a special feature of gluon 
jets which may not just be built from a long qq colour triplet string. 

While the study of spectroscopy did not provide a definitive evidence for a gluonic 
component in the spectrum, the inclusive production of the leading clusters in a gluon 
jet appears to be different from what is expected within a model which reproduces 
otherwise the main feature of the multiparticle final states in e^e~ annihilations. The 
lack of better signals from resonances may indicate that in the selection procedure for 
rapidity gaps particles from a resonance have been lost. Also the rapidity gap with 
Ay ^1.5 may be too small for significant results. 

9.2. Proposal for further measurements at the LHC 

There are some advantages to follow this line of studies at pp colliders, especially at the 
LHC [195]: 

• there is access to jets of much higher energies in comparison to LEP, therefore larger 
rapidity gaps are available with better separation of the leading cluster; 

• quark and gluon jets can be analyzed at comparable energies in the same 
experiment; 

• the gluon jet sample could be much larger than before. 

At pp colliders quark and gluon jets can be produced in 2 — )• 2 processes of leading 
order: quark jets in — )■ 7 + jet + X with subprocess qg — > 7g dominant at the lower 
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Figure 25. Mass density distributions for the leading cluster Mi^ading of charged and 
neutral particles; left: DELPHI results for gluon jets (JET3 of lowest energy in 3 jet 
events (sample 1) and in bbg events (sample 2)) and for quark jets (JETlof highest 
energy), all with charge Q leading — (figure from [200]); right: OPAL results for all 
mass combinations and for clusters with charges (H — ) and (H 1 — ), figure from [201]. 



Table 13. Purities (in%) of quark and gluon jets at Tevatron [204] and at LHC [194], 
taken from [195]. 

Pt xt g in di-jet q in 7+ jet 

Tevatron 1.8 TeV 50 0.056 60% 75 % 

LHC 7 TeV 200 0.057 60% 80 % 

50 0.014 75% 90 % 

800 0.229 25% 75 % 



Pt and gluon jets in di-jet events also at the smaller p^. The rates can be calculated 
from parton distribution functions and parton-parton cross sections, the purities depend 
on primary energy and p^, for some sample calculations see table 13. For quark jets a 
purity of 80% seems sufficient for the study of leading systems (quarks fragment harder 
than gluons), for gluons a higher purity of > 90% is desirable. This can be achieved by 
triggering on large total transverse energy and the selection of 3-jet events. Similar to 
the case of e^e'^ — )■ 3 jets the lowest momentum jet is most likely a gluon jet from QCD 
Bremsstrahlung {qqg, qgg or ggg). Estimates suggest that purities of 95% for gluon jets 
are within reach [195]. Studies at LHC could be important in two directions: 

9.2.1. Leading clusters with larger rapidity gaps. The energies of identified quark and 
gluon jets at the LHC could be larger by an order of magnitude as compared to LEP 
(typically < 25 GeV) and the rapidity gaps could extend up to Ay ~ 4. Then, the 
leading charges are expected closer to their asymptotic values Q leading = 0, ±1 with a 
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better sensitivity to an extra component with Q leading — 0- 

9.2.2. Identifying gluonic mesons through [qq]-reference x-spectra. The observed mass 
spectra of the leading cluster beyond the rapidity gap did not show prominent resonance 
structures. As the rapidity gap selection cuts into the angular distribution of the cluster 
decay particles the resonance signal may be reduced. In an alternative approach, the 
mass spectra of particle clusters M(7r7r), M{KK), M(47r) . . . are measured at a given 
total energy fraction x of the cluster and the x-distribution of resonances is determined. 
Such resonance x-spectra have been determined in the past [205] , but not separately for 
quark and gluon jets. 

The identification of gluonic mesons gh by the classification as "leading" and "non- 
leading" appears feasible by the comparison of the x-spectra D^^{x) with the spectra 
of reference [qq] mesons of similar mass D^'i^\x). In comparison with the "leading" 
spectrum like m — )■ tt"*" a gluonic meson should be suppressed 

< Dl^ for large x. (123) 

Similarly, a gluonic resonance foigb) should be suppressed in a quark jet as compared 
to a reference [qq] resonance and, vice versa, a [qq] meson should be suppressed in a 
gluon jet, as summarized in table 14 (for tetraquark mesons we assume here a steeper 
fragmentation as for [qq] mesons according to the counting rules [206]). As an example. 

Table 14. Mesons are leading or suppressed in fragmentation x-spectra depending on 
whether they carry the jet initiating parton {q or g) as a valence constituent; gluonic 
mesons are detected by comparison to reference qq mesons 



meson 




quark jet 


gluon jet 
triplet neutr. 


octet ncutr. 


M ■■ 


fo {ref : p,f2} 


leading 


suppressed 


suppressed 


gh: 


fo {ref : p,f2} 


suppressed 


suppressed 


leading 


[Aq]: 


a,/o(980) (?) 


suppressed 


suppressed 


suppressed 



for a gluonic scalar meson fo{gb) decaying to tttt (candidates /o(500), /o(980), /qISOO)) 
we can take p(770) and /2(1270) as reference [qq] states. Then, 

at large x : Df°^^^\x) < Dp/^{x); D^/^^'^ > Dp/^(x). (124) 

As found by the OPAL Collaboration [207] the x-distribution of /o(980) and /2(1270) 
almost coincide, then there is no evidence so far for any constituent structure beyond qq 
of /o(980) from this study; a similar measurement of the corresponding x spectra in gluon 
jets could reveal a gluonic admixture of /o(980) if there appeared a hard component at 
large x. In this way, also mixed gb — qq states could be recognized. 

Gluonic components could appear in the spectra of (7r7r)° from /o(500)/(t, /o(980), 
/o(1500); in (47r)° from /o(1370)(?), /o(1500) and in {KK)'^ from /o(980), /o(1500) and 
/o(1710). As an example, if /o(1500) has a large glueball component, it should fragment 
harder than p or /2 in a gluon jet and softer in a quark jet. 
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10. [qq] multiplets and glueballs from symmetry relations: 
first results from Xc decays 

In this section we investigate the exclusive production of hadronic states Hi,H2 in 
decays of heavy Quarkonia with C — +1 through the intermediate process 



in particular in decays of Xc or rjc {xb ot rjb). This decay could be viewed as fragmentation 
of gluons g ^ Hi,H2 with energy fraction z — 1 with sensitivity to glueball production 
assuming that the mass of Xc is high enough for the gg intermediate state to dominate. 
Here, we explore in particular the symmetry relations for the set of decay rates if all Hi 
and H2 are taken from [qq] multiplets Mi and M2 respectively. These relations follow 
from the quark states being produced as flavour singlet. 

A nice example is given by the results in table 4 on Xc decays where the properly 
weighted decay rates into pairs of pscudoscalars follow a simple cquipartition pattern. 
The same pattern should be followed by all pairs of mesons Hi, H2 taken from the same 
multiplets Mi, M2 each, i.e. for the branching ratios B into the members of two nonets 
after summing over charges 



with obvious notation. In this way, the classification of hadrons into multiplets can 
be tested. At the same time any additional gluonic component should become visible 
by a deviation from the symmetry relations. This will be discussed next for the scalar 
multiplet. 

10.1. Scalar multiplet according to route 2 

In this case, the scalars ao(980), /o(980), Kq{1430) and /o(1500) belong to the same 
multiplet. For three channels results on Xco decays have been obtained by BES2 
[208, 169] which we list in table 15 (taking the values from PDG [4]) and we correct for 
unseen channels. We take into account furthermore only corrections to the rates due to 
phase space, but no corrections from quark masses or additional dynamical form factors 
according to our findings on Xc decays above within the accuracy of around 10%. 

The normalized decay rate for /q/o is found slightly higher but compatible with 
the rescaled KqKq rate (1.64 ± 0.73 vs. 1.0), as expected for members of the same 
multiplet; irQ(1430) is generally accepted as qq state. An increased rate for decay into 
pairs of /o's could be indicative of an additional glueball component. The experimental 
uncertainties are still rather large for definitive conclusions. It is interesting to note the 
similarity of the two decay rates 



[QQ] ^gg^ H1H2, 



(125) 



B : [QQ] TTiTTs : K1K2 : K1K2 : 771772 : 77^77^ = 3:2:2:1:1 



(126) 



5(Xco ^ /o(980)/o(980)) = (2.68 ± 1.21) x 10-^ 
B(Xco 77'(958)77'(958)) = (2.02 ± 0.22) x lO'^ 



(127) 
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Table 15. Branching ratios B of Xco(3415) (0++) into pairs of scalars K*{U30)'^, fo = 
/o(980), /g = /o(1500) and baryons for comparison (data from PDG [4]) (i) in observed 
channels, (ii) corrected for unseen decays (isospin, also i?(/o(980) — ?• nn) — 0.75±0.12 
[169] and B(/o(1500) tttt&iKK from PDG); (iii) B/{cq) with charge weight c and 
momentum q after rescaling to the value for Kq{14:30)^ Kq{1430)^ , together with the 
resp. theoretical expectations for process gg as in table 4. 



weight c 


2 


/o/o 

1 


/o/o 
1 


PP 


AA 


(i) B X 10^ obs 

(ii) B X 10^ corr 

(iii) B/{cq) resc 


n qq+O-40 

u.»»_Q 29 
9 90+0.90 


(7r±7rT) 
0.67 ±0.21 

2.68 ± 1.21 
1.64 ±0.73 


(7r±i^=F) 

< 0.05 

< 5.0 

< 5.2 


0.223 ±0.013 


0.33 ±0.04 


expect (R=1.0) 


1. 


1. 


1. 







which again is in hne with the two states being parity partners as in route 2 with the 
same qq mixing; the similar decay rates in J/ip ^ {rj' , /o(980) + V) with V = {(f), u) have 
been noted ahcady in [61] (see section 4.2). The observation (127) for a gluonic decay 
indicates that not only the quark structure but also the gluonic components (if any) are 
similar in rj' and /o(980). A higher accuracy of the /o rates arc certainly desirable. 

A critical test for the existence of a scalar qq multiplet along route 2 would be the 
observation 

BiXco ^ ao(980)ao(980)) : B{xco ^ i^o* ( 1430) i?*( 1430)) =3:4 (128) 

if all charge states are counted, or 1:1 for single charge states. The rates for /o/o and /q/q 
could possibly deviate from the symmetry expectations because of gluonic components. 

The glueball should appear through a signal from Xco /o(500)/o(500) where 
/o(500)/(T is the broad object peaking at 1 GcV. The data in Xco 47r by BES [208] show 
indeed a considerable broad tttt background. Besides /o(980) there is some enhancement 
between 1200 and 1500 MeV which could get contributions from /2(1270) and /o(1500) 
interfering destructively with the background so as to yield a drop of the mass spectrum 
around 1500 MeV. Whether this background corresponds to the broad /o(500) could 
possibly be found out by observing the interferences with the narrow states /o(980) and 
/o(1500). 

In the decay of Xc2 [208] a strong signal of decay into pp shows up. In addition, 
there appear some bands around mTTTr ~ 500 MeV; this could be a signal from 
Xc2 — ^ /o(500)/o(500) where the mass spectrum is squeezed towards low masses because 
of the threshold factor in cms momentum p. This could be a signal from the decay 
into glueballs. 

If /o(500) has a high gluon component its pairwise production should be much 
favoured over the pairwise production of Kq (800) /urn Xco and Xc2 decays. 
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10.2. Scalar multiplet according to route 1 - tetraquark expectations 

In this case a, k, /o(980) and ao(980) are included in the lowest mass multiplet. If they 
are built from diquarks the decay would proceed according to the diagram in figure 3a 
but with production of diquarks instead of quarks (a contribution from figure 3b has 
been found negligible from the study of table 4). In this process at cms energy of 3.5 
GeV a formfactor for the diquark should play a role. This is actually observed for the 
production of baryon pairs which are produced from a primary diquark and a primary 
quark pair. These rates {pp, AA) are seen in table 15 to be an order of magnitude smaller 
than the rates for [qq] scalars like KqKq. One would therefore expect that scalars made 
up from two diquarks would be produced with a smaller rate yet. This appears to 
contradict the large rate of /o/o vs. KqK^ in table 15. 

Furthermore for route 1, one should observe the decays aa and kk with the 
same rate as /o/o and ooOq for particular charge states. In their analyses the BES 
Collaboration finds indeed a signal from kk of the order of K^Kq] they see signals with 
a but no signal from era where they take a according to the observation in J/ip ^ UTrn 
with a peak at low mass. The results on a and k are considered finally as too uncertain 
to be transformed into decay rates. In Xc2 decays there seems to be the channel aa 
as discussed above, but there is no evidence for the channel /o(980)/o(980). It will be 
interesting to see, whether the branching ratios follow 

B{xc ^ : B{xc ^ (7(7) = 4 : 1 (129) 

summed over all charge states, which is a result opposite to route 2 above. 

We conclude that the further study of Xc decays with higher statistics (or even of 
Xb decays) has a high potential to clarify the structure of the lowest scalar multiplet and 
the potential presence of gluonic mesons. The same, of course, is also true for multiplets 
with other quantum numbers, for example heavier pseudoscalars. 

11. Summary and conclusions 

In a summary of our results we try to answer some questions and we provide a minimal 
mixing scenario for the scalar mesons with glueball which takes into account some 
relevant findings. 

11.1. Is there any experimental evidence for the existence of glueballs? 

We have studied different scenarios for a scalar meson spectroscopy which includes the 
hghtest glueball of QCD. None of the schemes considered is without problem. It remains 
unsatisfactory, if the existence of the glueball as a fundamental prediction of QCD, at 
the end depends on whether such a controversial state like /o(1370) does exist. We have 
therefore looked after other criteria as well. 

There appears to be one observation which can be considered as a direct hint for 
the existence of glueballs. That is the observation of several LEP experiments that the 
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leading system in a gluon jet has an excess of neutral charge, in one measurement by 
40%, if compared with a standard Monte Carlo program. There is a new chance at LHC 
to follow this line of investigation further by obtaining a larger sample of gluon jets 
with higher energy to study jets with a larger rapidity gap. This should make the effect 
clearer and could possibly give a hint towards the relevant glueball mass. Future studies 
are in a "win- win" situation: either they provide better evidence for gluonic mesons, or 
else they show a deficiency of the "standard Monte Carlo's" with their hadronization 
models which would be important to understand as well, as these Monte Carlos are the 
backbone of the phenomenological studies with jets. 

11.2. Are there supernumerous states in the scalar sector? 
This question is related to the selection of scalar nonets. 

i) Extra state /o(1370). With the additional states /o(1500) and /o(1710) nearby there 
is one state too much (discussed as roTitc 1 and 3). Here our analysis has given a clear 
answer. We have not found convincing evidence for /o(1370) in any of the 2-body decay 
channels in the study of a large variety of reactions. In particular, elastic tttt scattering 
phase shift data put an upper limit of < 10% for its branching ratio into vrvr. There 
is room for further studies in the 47r channel, but so far the results are controversial. 
Further experimental investigations are desirable. 

ii) Extra state /o(500). In this scenario the i^*°(800)/fi; is not needed. There is some 
enhanced interaction at low Xtt energy but the Xtt phase shifts saturate below 40°, 
very unusual for a resonance. A possible clue could come from the study of symmetry 
relations for decay rates in Xc decays. 

11.3. Are there hints towards a gluonic component in isoscalar mesons? 

Our analysis of the production and decay pattern has shown that /o(1500) can hardly 
be described as a pure qq state, a particular solution with a gg contribution is suggested 
depending on the mixing angle. For /o(980) there is no hint for a strong gluonic 
component from the 77 decay width. The state /o(500) couples to strange and non- 
strange quarks but is not "flavour blind" so it could be a mixed gluonic-quarkonic state. 
A promising tool for further clariflcation of the flavour structure of isoscalars and their 
mixing angles are the study of Dg and especially of Bg decays. 

11.4. Future studies of symmetry relations in Xc decays 

Besides the study of leading clusters in gluon jets we find the study of decays of heavy 
quarkonia like Xco, Xc2 into pairs of scalar particles particularly interesting. This follows 
from the success of the symmetry relations for pairs of pseudoscalar particles. Such 
relations could establish which particles form a multiplet. Deviations from symmetry for 
the isoscalars would be a hint for glueball. The present results are consistent with /o(980) 
and ir*°(1430) to be in the same multiplet, /o(980) as tetraquark state is disfavoured. 
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Although the treatment of particles with a very large width is difficult, there may be 
some clues on a and k and their multiplets or gluonic component from such studies. 

11.5. A minimal scenario for scalar quarkonium-glueball spectroscopy 

Finally we present a solution which takes into account the main results on the flavour 
structure of the /o mesons discussed in this report. It is based on the qq nonet of route 2, 
i.e. it includes /o(980), ao(980), K*^{U30), /o(1500) with a flavour mixing similar to 
the mixing of ri,ri'. The near singlet flavour structure of /o(980) has been strongly 
supported with a mixing angle around 30° (as compared to 42° for the pseudoscalars) . 
Also the relative negative sign of the ss component within both /q states is established. 
In that case, /o(1500) needs a contribution from gg to explain the decay branching ratios 
(tttt, KK, Tjr}) with their small errors. At the same time, /o(500) does not look like 
a pure glue state, as it doesn't seem to realize flavour symmetry in the decays. The 
interpretation of 77 decays, sometimes with large errors, is not without controversy 
because of the model dependence. Our flndings on the flavour structure suggest as a 
minimal solution the existence of a glueball which is strongly mixed into /o(500 — 1000) 
and /o(1500) whereas /o(980) remains a quarkonium state. We recall that in the physical 
region the maximal strength of tttt scattering is near 1000 MeV (and not 500 MeV) with 
a comparable width. In this sense the broad /o(500 — 1000) and the narrow /o(1500) 
overlap in mass, a situation reminiscent of the p — uu system, apart for the very different 
widths. This suggests a mixing 



We prefer for /o(1500) the solution tan^c > and especially with our value of 4>, 



This may be compared with maximal qq — gg mixing at (f)a ~ 45° while the model [61] 
corresponds to (pc = 0. The result on (f)G assumes a relation between the decay couplings 
{go ~ Gq); more generally, one may introduce an extra parameter Tgq = Go/go, for now 
rco ~ 1- The flavour structure of /o(1500) with dominant gg and qq octet components 
is also comparable to the models [35, 51] (see (26)) but with a very different multiplet 
structure. 

Specific experimental tests of the presented scenario include the better 
determination of mixing angles of the /o mesons and their gluonic components. Also the 
amplitude signs lead to observable effects. A special property of /o(500 — 1000) in this 
scheme is the negative coupling of the ss component with respect to the one of /o(980) 
which is indicated in the Argand diagram of figure 13. 

A new possibility for glueball spectroscopy is opened up by the gluonic coupling of 
/o(500)/(T and its dominant decay into nn. Because of this coupling a higher mass 



1/0(500 - 1000)) = sin0a|gg) - cos ^dgg) 
|/o(1500)) = cos (j)G\qq) +sm(t)G\gg) 

\qq) = cos (l)sc\nn) — s\n(f)sc\ss) 



(130) 
(131) 
(132) 



(t>sc = 



(30 ± 3)°, ~ 35°. 



(133) 
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glueball should have a favoured decay into aa. In case of /o(1500), which has a 
probabihty 1/3 for a gluonic component according to the above simple estimate, there 
is a decay BR of around 13% into a a according to table 12; the evidence is particularly 
clear and unique for the decay channel o — )> 7r°7r°. 

As an interesting aspect of the scenario presented for scalar spectroscopy we note 
the similarity to the recent solution from QCD sum rules for the scalar sector which 
suggests two gluonic contributions at 1000 and 1500 MeV, but with strong correlation. 
Also there is one group with a similar result from lattice calculations, but alternative 
views exist as well. 
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